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Volume  V  deals  with  tlie  following  topics; 

1)  detectlou  in  an  snisotronlc  noise  field 

Earlier  studies  of  passive  detection  in  an  anisotropic  noise 
environment  were  only  concerned  with  the  anisotropy  caused  by  a 
single  plane  wave  interference  (Volumes  III  and  TV) .  The  present 
volume  presents  analyses  of  detection  in  a  noise  field  dominated 
by  several  plane  wave  interferences  or  by  a  single  spatially 
distributed  interference.  Conventional  as  well  as  optimal 
detectors  are  consiaered. 

2)  Passive  tracker  accuracy 

a)  One  study  examines  the  effect  of  a  plane  wave  interference 
on  the  performance  of  a  split  beam  tracker.  Conventional  as 
well  as  null  steering  types  of  beam-formers  are  ar  iyzed. 

The  contribution  of  the  interference  to  the  measured  bearing 
error  is  decomposed  into  a  systematic  and  a  random  component. 
Factors  which  affect  the  relative  magnitude  of  these  components 
are  discussed. 

b)  A  second  study  initiates  an  effort  to  set  absolute  lower 
bounds  on  the  bearing  accuracy  attainable  with  a  given  array 
in  a  given  noise  environment.  Only  the  simplest  case  (two 
element  array.  Independent  noise)  is  discussed  here. 

3)  Active  receivers  using  replica  correlation 

The  performance  of  a  simple  replica  correlator  is  compared  with 
that  of  a  similar  instrumentation  using  clipped  (binary)  hydrophone 
Detc  tion  as  well  as  range  and  Doppler  measurement  are 


iii 


outputs . 


examined .  One  finds  that  large  clipping  losses  can  occur  when 
the  target  is  moving  rapidly.  In  most  other  situations  the  clipping 
loss  is  small. 

4)  Adaptive  Signal  Processing 

In  situations  where  processor  design  is  hampered  by  lack  of  adequate 
knowledge  concerning  signal  or  noise  statistics  one  can  employ 
stochastic  approximation  techniques  to  cause  the  processor  to  approach 
an  optimum  configuration.  This  procedure  is  used  to  adjust  a  tapped 
delay  line  filter  for  operation  in  a  noise  environment  with  unknown 
spectral  properties.  Conditions  of  convergence  and  rates  of  conver¬ 


gence  are  examined. 
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FOliEWORD 


This  report  is  the  fifth  in  a  series  desciibing  work  performed  by  Yale  University 
under  subcontract  to  Electric  Boat  division  of  General  Dynamics  Corporation.  The 
report  covers  the  period  from  1  July  IfXhi  to  1  July  1967.  An  imclassified  supple¬ 
ment  to  this  volume  has  been  bound  as  a  separate  document  (U417-68-079).  Electric 
Boat  is  prime  contractor  of  the  SUBIC  (Submarine  Integrated  Control)  Program 
under  Office  of  Naval  Research  contract  NOnr  2512(00U  LCDR.  E.W,  Lull,  USN,  is 
Project  Officer  for  ONR;  J  W.  Herring  is  Project  Manager  for  Electric  Boat  division 
under  the  direction  of  Dr.  A.  J.  vanV.^oerkom. 


I. 


Introduction 


The  following  is  a  summary  of  work  performed  under  contract 
8050-31-55001  between  Yale  University  and  the  Electric  Boat  Company  during 
the  period,  1  July  1966  to  30  September  1967.  More  detai.led  discussions 
of  the  results  as  well  as  their  derivations  are  -.ontalned  in  a  series  of 
six  progress  reports  which  are  appended. 

Several  studies  reported  in  earlier  volumes  of  this  series  were 
concerned  with  passive  detection  in  the  presence  of  strong  interference  from 
a  point  source.  Tlio  general  subject  of  passive  sonars  operating  in  an 
anisotropic  noise  environment  is  pursued  further  in  the  present  volume. 

The  effort  reported  here  has  taken  two  new  directions: 

1)  Anisotropies  are  no  longer  attributed  to  a  single  interfering  plane 
wave.  Environments  containing  several  point  sources  of  interferenct 
or  a  spatially  distributed  interference  are  studied.  Conventional 
as  well  as  optimal  detectors  are  analyzed. 

2)  The  effect  of  a  single  plsne  wave  interference  on  tracking  accuracy 
is  examined. 

The  problem  of  tracking  accuracy  Is  also  considered  in  another  context. 
An  attempt  is  made  to  divorce  the  accuracy  problem  from  particular  instru¬ 
mentations  and  to  set  absolute  bounds  on  the  bearing  accuracy  att.alnable  by 
processing  the  outputs  of  a  set  of  hydrophones  operating  in  a  specified  nolst 
environment.  Only  the  simplest  rossible  case  (two  hydrophones,  noise 
Independent  from  phone  to  piione)  is  presented  here. 

In  addition,  th'.s  voliirae  continues  the  study  of  active  sonar  systems 

is 

Initiated  in  Volun.e  IV.  The  specific  problem  considered  here/the  effect  of 
clipping  on  the  pertormance  of  replica  correlators.  Results  are  obtained 
for  a  wide  class  of  signal  waveshapes  and  for  environments  dominated  either 
by  reverberation  or  by  ambient  noise. 


Finally,  an  effort  is  initiated  to  deal  with  the  signal  detection  and 
extraction  problem  in  a  noise  environment  whose  statistical  properties  are 
largely  or  wholly  unknown.  This  leads  to  the  study  of  adaptive  processing 
procedures.  Only  preliminary  results,  based  on  the  method  of  stochastic 
approximation,  are  presented  in  this  volume. 


II.  Detection  In  an  Anls o t roplc  H olse  Environment 

Report  No.  30  deals  with  the  performance  of  a  conventional  detector 

operating  in  a  noise  environment  dominated  by  several  point  sources  of 

Interference  or,  as  a  limiting  case,  by  an  interference  source  spatially 

distributed  over  some  finite  angle.  The  results  indicate  that  the 

performance  degradation  due  to  such  complex  Interference  sources  is  less 

serious  than  that  caused  by  a  single  point  source  of  interference  with  the 

same  total  power.  For  a  linear  array  of  M  equally  spaced  hydrophones 

the  maximum  differential  amounts  to  lO.loq  v^K/3  db  of  equivalent  input 

in 

signal-to-noise  ratio,  .■'ais  is  precisely  the  performance  differential 

between  .1  conventional  array  operating  in  a  noise  environment  Independent 

from  hydrophone  to  hydrophone  and  a  similar  array  operating  in  a  noise 

environment  of  equal  j/ower  but  criginating  largely  from  a  point  source  of 

interference.  One  therefore  suspects  that  a  smooth  transition  will  take 

place  from  the  cast  of  a  single  point  source  of  interference  to  that  of 

isotropic  noise  as  the  number  of  interferences  increases  and  their  locations 

become  more  uniformly  distributed  in  space.  Numerical  computations 

generally  conflnn  this  Inference.  For  several  closely  spaced  Interferences 

all  relatively  remote  from  the  target  in  bearing  one  finds,  not  surprisingly, 

that  the  configuration  is  equivalent  to  a  single  interference.  As  the 

spacing  between  interferences  increase?  the  performance  index  rises  quickly 

to  Its  asymptotic  value.  For  the  situat  ons  considered  in  the  computations 

the  asymptotic  improvtroent  Is  somtwhaU  less  than  the  figure  of  lO-log  /^/3  db 

10 

quoted  above  One  major  reason  for  this  discrepancy  is  the  fact  that  the 
postulated  spectra  (falioEf  with  second  power  ni  frequency  above  5000  cps) 
are  not  sufficiently  broad  to  yield  noise  independent  from  hydrophone  to 
hydrophone,  in  the  isotropic  limit. 


T 


Qualitatively  similar  results  ere  obtained  for  the  case  of  a  distributed 
interference.  Again  one  finds  that  the  performance  index  quickly  approaches 
its  asymptotic  maximum  as  the  interference  spread  grows.  This  approach 
becomes  more  rapid  as  the  number  of  hydrophones  increases  (and  hence  the 
beam-v.’ldth  decreases,  since  fixed  spacing  between  hydrophones  is  assumed). 

Report  No.  31  studies  a  likelihood  ratio  detector  operating  in  an 
anisotropic  environment  similar  to  the  one  Jest  discussed.  It  had  been  shown 
in  earlier  work  (Volume  III)  that  interference  from  a  single  plane  wave 
resulted  in  a  performance  degradation  noi  exceeding  the  degradation  caused  by 
the  loss  of  one  hydrophone  in  the  absence  of  the  interference.  A  similar 
statement  remains  true  in  the  presence  of  several  Interferences,  i.e.,  the 
effect  of  R  interferences  can  be  eliminated  at  a  cost  not  exceeding  R 
hydrophones.  However  .  it  appears  that  this  bound  on  performance  degradation 
is  often  quite  pessimistic.  Thus  interferences  in  close  angular  proximity 
of  each  other  have  the  effect  of  a  single  interference  and  can  be  dealt 
with  at  a  sacrifice  of  only  one  hydrophone.  Only  when  all  interferences 
are  widely  separated  from  each  other  and  from  the  target  can  the  loss  figure 
approach  R  hydrophones  (and  then  only  in  a  strongly  interference  dominated 
environment) .  A  spatially  distrlbuTed  interference  may  again  be  interpreted 
as  the  limiting  version  of  a  large  number  of  closely  spaced  individual 
Interferences.  Those  may  then  be  grouped  into  clusters  more  or  less 
equivalent  to  single  point  sources  of  interference  so  that  serious  losses 
occur  only  if  the  number  of  clusters,  i.e.,  the  total  angular  spread  of  the 
interference,  U-  large.  Since  analytical  evaluation  of  the  likelihood  ratio 
dutector  leads  to  expressions  that  are  formally  simple  but  practically 
difficult  to  evaluate  and  intirprct,  many  of  the  conclusions  are  based  on 
numerical  computations.  Arrays  of  different  geometries  (circular  and  linear) 
'.nd  a  variety  cf  interferen  e  configurations  are  considered. 
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III.  Tracking  Accuracy 


Report  No.  29  analyzes  the  performance  of  a  split  beam  tracker  operating 
in  an  environment  consisting  of  ambient  noise  (assumed  independent  from 
hydrophone  to  hydrophone)  and  a  plane  wave  interference.  Signal,  noise  and 
interference  are  assumed  to  be  stationary  Gaussian  proC’-sses  with  similar 
power  spectra.  Twc  types  of  beam-formers  are  considered: 

1)  Convcntlnr.al  beam-formers  (in  each  array  half:  delay  for  alignment 
with  target,  then  add). 

2)  Null  steering  (in  each  array  half:  steer  on  interference,  subtract 
hydrophone  outpr.ts  pairwise,  then  beam-form  on  target  and  add) . 

In  each  case  one  cf  the  btams  is  shifted  90°  in  phase  relative  to  the  other, 
after  which  multiplication  and  smoothing  completes  the  tracking  procedure. 

In  the  conventional  tracker  plane  wave  interferences  (or  other  spatial 
asymmetries  in  the  noise  field)  contribute  to  the  tracking  error  through  two 
distinct  mechanisms. 

a)  The  null  of  the  average  tracker  output  is  shifted  away  from  its 
nominal  location  by  an  amount  depending  on  the  proximity,  strength 
and  spectral  properties  of  the  Interference  (systematic  error) . 

b)  The  output  fluctuation  is  increased  by  an  amount  depending  primarily 
(except  for  interferences  close  to  the  target  in  bearing)  on  the 
interference  power  (random  error). 

The  null  steering  tracker  eliminates  the  complete  time  function  of 
interference  (at  least  ideally)  and  therefore  removes  both  sources  of  error. 

Whether  tracking  problems  caused  by  noise  field  asymmetry  warrant  the 
use  of  instrumentations  more  complex  thau  the  conventional  split  beam 
tracker  clearly  depends  on  the  total  additional  error  [  a)  plus  b)  ].  However, 
the  choice  of  remedial  measures  (8pectr\im  shaping,  null  steering,  more 


powerful  adaptive  procedures,  etc.)  might  depend  to  a  considerable  extent 
on  the  relative  magnitude  of  c-rrore  a)  and  b)  .  The  following 
considerations  appear  relevant: 

1)  Only  when  the  noise  field  Is  interference  lomincted  does  the 
interference  contribute  significantly  to  the  random  error.  (With 
a  linear  a-rray  of  2M  equally  spaced  hydrophones  the  environment 
is  interference  doiainatea  when  the  interference  to  amolent  noise 
ratio  exceeds  /(4/3)M  ). 

2)  Even  when  the  noise  field  is  interfcrevice  dominated  the  systematic 
i.  rror  tends  to  exceed  the  random  error  in  many  practically 
Interesting  situations. 

3)  A  certain  amount  of  control  can  be  exurted  over  the  systematic 
error  by  shaping  the  spectrum  of  eacli  ch.annel  prior  to 
multiplication . 

If  Che  primary  source  of  difficulty  is  systematic  error,  only 
st,.tistical  properties  of  the  asyminetrical  noise  component  are  required  for 
correction.  These  may  not  be  known  a  priori  so  that  measurement  and  at 
ler.st  a  primitive  form  of  adaptation  may  be  indicated.  However,  If  one 
wishes  to  deal  with  the  random  component  of  error,  one-  must  obtain  an 
estimate  of  the  actual  Intcrf e>rence  time  function  and  the  need  for  adaptation 
becomes  more  fundamental.  Even  in  ch.-it  case  one  should,  of  course,  still 
use  all  available  a  priori  information.  Thus,  if  one  knows  chat  the  noise 
field  asymmetry  is  caused  by  an  interfering  plane  wave,  one  need  only 
tieasure  the  interference  bearing.  Then  the  relatively  simple  null  steering 
procedure  solves  the  problem.  In  the  absence  of  such  strong  information 
concerning  the  cause  of  noise  field  asymmetry  one  nay  be  forced  into  more 
elaborate  adaptive  tecliniques. 


Report  No.  32  considers  the  tracking  probli.r  from  a  much  more  general 
point  of  view.  The  basic  aim  is  to  deterr.ilrie  the  best  passive  bearing 
accuracy  attainable  with  a  given  array  in  a  given  noise  environment,  without 
making  any  prior  assumptions  about  data  processing  procedures.  Report  No.  32 
deals  only  with  the  simplest  possible  case,  that  of  a  two  element  array  with 
noise  independent  f~om  hydrophone  to  hydrophone.  In  that  situation  a  simple 
crosscorrelator  attains  the  Cracicr-Rao  lower  bound  of  bearing  error  and  is 
therefore  an  optimal  bearing  estimator.  A  two  element  split  beam  tracker 
using  differentiation  to  achieve  the  90°  phase  shift  betvieen  channels  is 
equivalent  in  performance  to  the  autocorrelator.  If  a  pure  90°  phase  shift 
is  used  in  place  of  the  differentiation  there  is  a  slight  degradation  of 
perforaance.  The  rms  bearing  error  of  the  optimal  instrumentations  varies 
as  (S/N)  ^  for  S/N  <<  1  [S/N  =  input  slgnal-to-noise  ratio)  and  as 
(S/N)”^  for  S/N  >>  1  . 
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IV.  Clipped  Ri-plice  Com-lcitorg 


Report  No.  31  u.'cnr.lriv.s  the  performance  of  an  .active  sonar  receiver 
using  replica  corre'^tlon.  Detection  r.s  well  as  range  ind  Doppler  neasure- 
nent  are  considertc.  Th^.  output  of  each  array  decent  is  clipped  prior  to 
conventional  bcanforclng  and  the  subject  of  primary  interest  is  the 
performance  degrnd.'-.tlon  due  to  clipping.  Ambient  noi.'Je  as  well  as  reverbera¬ 
tion  limited  environments  are  studied,  with  primary  emphasis  on  the  latter. 
The  rcverbe‘ration  model  is  the  one  developed  in  Volume  IV  (stationary, 
independent,  Poisson  distributed  point  scatterers)  . 

The  general  conclusions  m.ay  be  summarized  as  follows:  Serious  clipping 
losses  can  arise  In  a  reverberation  United  envi*-onirient  when  the  target  is 
moving  rapidly  enough  to  shift  the  target  return  slmost  entirely  out  of  the 
reverberation  band.  In  such  cases  the  undipped  detector  h.-.s  an  output 
essentially  free  from  reverbe.r'ition  in  the  signal  band,  whereas  the  clipping 
opeiatlon  transfers  part  of  the  reverberation  power  into  the  signal  band. 

In  other  words,  clipping  destroys  much  of  the  potential  slgnal-to-noise 
advantage  of  a  rapidly  coving  t.arget.  With  this  exception  the  clipping  loss 
appe.acs  to  be  quite  small  in  most  practically  interesting  situations.  If 
one  defines  the  clipping  loss  R  as  the  ratio  of  the  output  slgnal-to-noise 
ratios  with  .and  without  clipping,  one  can  show  under  fairly  general 
conditions  that  R  >  0.89  (equivalent  to  a  loss  of  about  1  db  of  input 
signal- to-noisc  ratio).  The  primary  requirement  for  this  statement  to  be 
true  is  that  the  tr.'nsmit  ted  signal  be  n.arrowb.ind  in  scu.ic  meaningful  sense 
(in  <a  reverberation  United  environment  It  is  sufficient  -  but  by  no  weans 
necessary  -  that  the  wavelength  of  the  highest  modulation  frequency  be  l.arge 
compared  with  the  array  dimensions).  In  the  absence  of  such  minimal 
requirements  one  can  construct  pathological  examples  which  yield  values  of 
R  arbitrarily  close  to  zero. 


The  structure  and  performance  of  the  optimum  detectors  discussed  in 
earlier  volumes  of  this  series  depends  critically  on  prior  knowledge  of 
signal  and  noise  statistics.  It  has  been  pointed  out  that  such  knowledge  la 
likely  to  be  Incomplete  at  beat.  Several  reports  have  postulated  that  certain 
noise  parameters  (.e .  g .  total  noise  power)  were  unknown  and  have  studied  the 
resulting  detection  problem,  arriving  at  primitive  forms  of  adaptation  to 
the  noise  field. 

Report  No.  34  represents  the  first  attempt  in  this  series  to  deal  with 
a  truly  adaptive  situation,  one  in  which  there  I s  an  absolute  minimum  of 
available  a  priori  information.  The  report  is  of  an  introductory  nature. 

It  seeks  to  define  the  problems  and  create  the  mathematical  framework  for 
later  studies.  Specific  results  ore  obtained  only  for  a  signal  processor 
using  a  single  hydrophone,  but  the  procedures  discussed  generalize  without 
difficulty  to  array  processing  problems. 

In  the  situation  analized  here  the  output  of  the  single  hydropnone 
is  passed  through  a  tapped  delay  line  filter  whose  tap  weights  are  adaptively 
adjusted  to  minimize  the  mean  square  error  between  the  filter  output  and 
the  signal  component  of  the  input.  Stochastic  approximation  is  used  as 
the  basic  technique  of  adjustment. 

It  is  clear  that  detection  is  impossible  in  the  total  absence  of  any 
information  concerning  signal  or  noise.  Three  cases  short  of  such  total 
ignorance  are  considered;  1)  The  signal  waveshape  is  known,  the  noise  is  a 
stationary  stochastic  process  with  unknovm  statistical  properties  2)  Signal 
and  noise  are  stationary  stochastic  processes,  the  signal  spectrum  is 
known,  the  noise  spectrum  is  not  3)  Signal  and  noise  are  stationary  stochastic 
processes,  the  noise  spectrum  is  knovm,  the  signal  spectrum  Is  not. 


In  each  case  one  finds  that  proper  use  of  the  stochastic  approximation 
technique  generates  a  filter  whicn  converges  to  the  optimum  Wiener  filter. 
Conditions  of  convergence  and  rates  of  convergence  are  similar  for  ’’he 
three  cases.  A  convenient  choice  of  the  gain  parameter  in  the  stochastic 
approximation  algorithm  causes  the  mean  square  error  to  converge  to  its 
minimum  with  the  approximate  first  power  of  time.  Certain  preliminary  re¬ 
sults  dealing  with  the  application  of  stochastic  approximation  techniques 
to  non-stationary  noise  fields  are  also  presented. 
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Summary 


The  report  analyzes  the  effect  of  a  plane  wave  interference  on  the 
performance  of  a  split  beam  tracker.  The  performance  of  a  conventional 
tracker  is  compaied  with  that  of  a  tracker  designed  to  null  the  interfer¬ 
ence  prior  to  beam  forming  on  the  target.  The  receiving  array  is  assumed 
to  be  linear  with  equally  spaced  hydrophones.  Signal,  ambient  noise  and 
interference  are  assumed  tu  be  statistically  independent  Gaussian  random 
processes  with  power  spectra  of  the  same  form.  The  observation  time  T 
is  assumed  to  be  large  compared  with  the  correlation  time  of  signal, 
ambient  noise  and  interference,  but  short  enough  so  that  target  and  inter- 
fi'rence  bearings  do  not  change  significantly  in  T  seconds.  For  computa¬ 
tional  simplicity  the  ambient  noise  is  assumed  to  be  statistically  inde¬ 
pendent  f ror,  hydrophc.te  to  hydrophone.  The  following  results  are  obtained 

1)  Ihe  track!’ g  error  consists,  in  general,  of  two  parts: 

■i)  Systematic  error  is  measured  by  the  displacement  of 

the  target  null  of  the  average  cracker  output  from  the  true 
target  bearing.  It  is  due  to  asymmetry  in  the  noise  field, 
caused  here  by  the  interference, 
b)  Random  error  is  the  fluctuation  of  the  tracker  output  about 
its  average  value.  It  is  due  to  the  finite  smoothing  time  T 

2)  The  nulling  tracker  completely  eliminates  the  interferei.ee.  With 
the  spatially  symmetrical  ambient  noise  postulated  in  the  analysis, 
its  output  is  therefore  free  from  systematic  error  and  exhibits  a 
random  error  depending  only  on  the  ambient  noise. 


A-i 


3) 


4) 


3) 
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Thi-  convention,"!  tracker  has  both  systematic  and  random  error 
components.  The  ratio  of  systematic  to  random  error  is  greatest 
for  spectra  whose  Fourii r  tr'nsforms  decay  slowly  (e.g.,  white 
spectra)  ai  smallest  for  spectra  with  rapidly  decaying  Fourier 
transforms.  However,  even  in  the  latter  case  the  systematic  error 
tends  to  exceed  the  random  error  (often  by  a  large  factor)  under 
most  reasonable  operating  conditions,  as  long  as  the  environment 
is  interference  dominated.  For  interferences  well  separated  from 
the  target  the  borderline  between  interference  dominated  and 
ambient  noise  dominated  operation  is  reached  at  an  ambient  noise 
to  Interference  ratio  of  / (4/3)M  .  M  is  the  ntimber  of  hydro¬ 
phones  in  each  half  of  the  array. 

At  low  interference  to  signal  ratios  the  systematic  error  of  the 
conventional  tracker  rises  linearly  with  the  interference  to 
signal  ratio.  The  slope  of  this  rise  depends  strongly  on  the 
spectral  shape,  being  largt;st  for  spectra  with  slowly  decaying 
Fourier  transforms.  Spectrum  shaping  can  be  accomplished  by 
insertion  of  appropriate  filters  into  each  channel  of  ttie  tracker. 
When  tlie  interference  to  signal  ratio  increases  beyond  a  certain 
point  the  systematic  error  increases  rapidly.  Soon  thereafter 
the  target  null  disappears  entirely  and  tracking  becomes  impossible. 
This  may  happen  at  interference  levels  at  which  the  random  error 
is  not  at  all  excessive.  Even  at  interference  to  signal  r.atios  as 
low  as  5  the  systematic  error  can  easily  amount  to  a  major  fraction 
of  a  degree. 

With  fixed  hydrophone  spacing  the  rms  random  error  of  the  conven- 

T  3 / ?  3/2 

tional  tracker  varies  as  T  ut  '  M  I  in  an  interference 


CONFiDENTiAL 


L)(;  ,\ol  idi’Y 


dominatlid  environment  and  as  T  ^  u  M  ^  N  in  an  ambient 

o 

noise  dominafac  environment,  cj  is  the  bandwidth  of  the  spectrum, 

o 

I  the  average  interference  power  and  K  the  average  noise  power. 

For  fixed  array  length  (hydrophone  spacing  inversely  proportional 

_L 

to  M  )  the  M  dependence  is  M  (int._rference  dominant)  and 
M  ^  (ambient  nols  dominant)  respectively. 

6)  The  ms  error  of  the  nulling  tracker  is  approximately  the  same  as 
that  of  the  conventional  tracker  in  on  ambient  noise  dominated 
environment.  In  an  interference  dominated  environment  the  ms 

e. '•or  'if  the  conventional  tracker  is  larger  by  the  (actor 
/8/¥  /m  (I/N)  . 

7)  It  is  apparent  from  2)  and  6)  that  one  might  employ  interfer¬ 
ence  nuLiIng  for  two  reasons 

'.)  To  eliminate  systematic  error 
b)  To  reduce  random  error 

Since  nulling  achieves  significant  random  error  reduction  only  in 
a  strorgly  interference  dominated  environment,  one  must  inquire 
whether  such  an  environment  occurs  suf I'lclertly  often  to  justify 
the  added  complexity  of  instrumentation,  if  it  does  not,  one  must 
further  inciuirc  whether  one  cannot  eliminate  system.atic  error  by 
procedures  much  simpler  to  implement  than  nulling.  This  appears 
distinctly  possible  because  sys ter...atic  error  depends  only  on 
average  parameters  of  the  interference  (bearing,  power,  etc.),  not 


on  details  of  the  interfere.-ice-  time  function. 
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1.  Introduction 

This  report  is  concerned  with  the  effect  of  interference  from  a  point 
source  on  the  pcrfonnance  of  split  beam  trackers.  Instrumentations  with 
and  without  provisions  for  null  steering  on  the  interference  are  analy^.ed. 
Because  interest  centers  on  the  Intekfcrcnce  problem,  the  simplest  possible 
assumptions  are  made  concerning  all  other  aspects  of  the  system.  Thus 
interference,  ambient  noise,  and  signal  are  assumed  to  be  independent 
Gaussian  random  processes  with  spectra  identical  over  the  processed  frequency 
band.  The  ambient  noise  is  regarded  as  statistically  i.  dependent  from 
hydr( ohone  to  hydrophone.  The  processing  array  is  linear  ar.d  "onsi'ts  of 
equally  spaced  elements. 

Varii  ^  possible  instrumentations  differ  considerably  in  their  detailed 

characteristics,  but  each  exhibits  an  average  output  z  which  varies  with 

ste..  ;ing  angle  6  roughly  in  the  manner  outlined  in  Figure  1,  at  least  for 

values  of  0  close  to  tic-  target  bearing  and  reasonably  remote  from  the 

interference  bearing  In  the  absence  of  .systematic  error,  z  passes 

chroiM',’..  zero  at  0  =  6  The  Instantaneous  tracker  output  is  the  value  of 

e  near  at  which  the  random  variable  z(f0)  (whose  mean  value  is  z) 

passes  through  zero.  Tiie  rttis  trackin,;  error  is  therefore  the  standard 

I  _ 

I  z 
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deviation  of  this  ztro  location.  As  long  as  this  standard  deviation  is 

small  compared  with  9^^  there  is  no  serious  danger  that  the  target  will  be 

lost.  On  the  other  hand,  if  the  rms  error  becomes  comparable  to  6,, 

M 

(i.e.,  if  2(t)  can  exceed  with  a  probability  that  is  not  negligible), 
sustained  tracking  is  no  longer  practical. 

Formal  computation  of  the  zero  distribution  of  a  random  process  such 
as  z(t:e)  (considered  as  a  function  of  0)  is  an  extremely  difficult 
problem.  However,  a  simple  approximation  can  be  obtained  for  the  situation 
of  primary  interest  here,  the  case  of  trackers  employing  smoothing  times 
large  compared  with  the  correlation  time  of  the  tracker  input  .  The 
assumption  of  large  smoothing  times  has  two  important  consequences: 

1)  The  tracker  output  z(t:0),  considered  ns  a  function  of  t,  is 

an  approximately  Gaussian  random  process.  Furthermore,  z(tj  :8j^)  ,z(t^ -62)  ) 
may  be  regarded  as  Gaussian  in  two  dimensions. 

2)  Smoothing  over  long  period  of  time  results  in  relatively  small 

scattering  in  both  amplitude  and  slope  of  z(t  .•6)--corisidered  as  a  function 
of  0 — about  the  averages  specified  by  z.  Sensitivity  considerations 
require  the  slope  to  be  reasonablv  large  near  6=0  ,  But  if  r-' 

d  w  '  ■  0  00 

is  a  sizable  positive  number  (as  suggested  by  Figure  1)  and  the  scattering 

By 

in  slope  is  small,  the  probability  of  — ■  assuming  a  ncgativi:  value  near 

0  U 

e  is  extremely  small.  In  other  words,  with  a  probaollity  close  to  unity, 
o 

z(t;6) — considered  as  a  function  of  6 — has  one  and  only  one  zero  in  the 
neighborhood  of  6^. 

^For  a  detailed  statement  of  the  required  conditions,  see  M.  Rosenblatt, 
"Some  Comments  on  Narrow  Band-Pass  Filters,"  Quarterly  of  Applied  Math., 

15 ,  No.  A,  January  1961. 
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As  soon  as  the  possibility  of  multiple  zeros  in  the  0  range  of 
interest  can  be  ruled  out,  it  becomes  a  relatively  simple  matter  tc 
calculate  the  zero  distribution  by  computing  the  closely  related 
probability  that  z(t-,6)  has  a  zerc  in  (6,6+d0).  Designating  the;  latter 
probability  by  P(0)  d0  one  obtains 


P(0)  d0  -  Pr{z(t;e+d0)  >  0}  -  Pr{z(f0)  >  0) 


(1) 


Hence 


P(0) 


30 


Pr{z(t;0)  >  0} 


J 


(2) 


In  view  of  the  Gaussian  nature  of  z(t;0) 


(z-z)‘ 
2a  ^ 


Pr(z(t,0)  >  0}  = 


/irt  a 


dz 


(3) 


where  o  is  the  standard  deviation  of  z.  Doth  z  and  a  are,  in  general, 
z  z  ’ 

functions  of  0.  Substituting  into  (2)  one  obtains 


P(e)  = 


fz-z)^ 

2o  ^ 
z 


30 


/2i 


dz 


/lit  a 


f 

/  dz  e 


0 


—  /•  2  3o 

z  >  z  -  z  3z  (z  -  z)  _ z 

G  ^  39  0  ^  30 

z  z 


(^) 


Straightforward  evaluation  of  the  inti,grr.l  leads  to  the  result 

-2 


P(0)  = 


1  ^ 

3z  z 

1 

1  +  erf 

7 

Bo 

z 

T=  e 

/  2tto 

z 

^30 

30  j 

2o 

z 

30 

(5) 
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where 


f  -  2  -y  A 

erf  X  : - /  c  ■'  dy 


If  the  standard  deviation  of  z  is  constant  over  the  6  range  of  Interest, 
3o 

“  0  and  Eq.  (5)  reduces  to 


_  1  z  37 

p(e)  =» -  e  — 

/2v  a  39 

z 


Sustained  tracking  is  clearly  feasible  only  if  the  typical  fluctuation  in 


indicated  0  is  small  compared  with  J -  6^]  (sec  Figure  1).  In  that 
case  0^  must  be  small  compared  with  z..  and  one  can  approximate 


z  =  +  —  (6  -  6  ) 

o  30  o 


z  and  ~  are  the  values  of  z  and  its  derivative  at  8  =  0  . 
o  30  0 

0 


Substituting  Into  Eq .  (7)  one  then  obtains 


P(0)  =  - ^ - expJ-  ~ 


®o  371 


V. 


Thus  the  zero  location  (the  indicated  bearing)  has  an  approximately 


z 

normal  distribution  with  mean  6  -  — — = —  and  standard  deviation 

o 

3^lo 

The  latter  is,  of  course,  the  ratio  of  fluctuation  to  sensitivity  at  the 

null,  the  figure  of  merit  frequently  used  to  measure  the  accuracy  of  null 
z 


iZ 

36  0 


seeking  devices. 


is  a  systematic  error  in  Indicated  bearing  due  to 


36! 


0 


such  factors  as  asymmetry  in  the  noise  field  or  imperfections  in  thc 
instrumentation . 

If  the  standard  deviation  of  z  varies  with  0  one  must  revert  to 

Iq  (5).  However,  for  most  cases  of  interest  in  the  present  study  it  will  be 

3o 

2 

found  that  the  contribution  to  Eq .  (5)  of  terms  involving  is  small,  so 

that  Eq .  (9)  givus  a  reasonable  approximation.  This  is  almost  evident  by 

inspection  of  Eq .  (5),  f.  r  if  one  multiplies  thi.s  oauation  by  36  it  read.= 

-2 


P(6)  36 


■/2ii 


1  2°z'  ' 

--  -  e 


az - 3o  - 


2a 


1  +  erf  — -  — 
/?  0 


3o 


(5a) 


Consider  valurs  of  z  satisfying  jzj  <_  o  ,  the  range  in  which  the 


.xponential  function  has  a  value  significantly  different  from  zero. 


mission  of  the  term 


smaller  than  the  change 


it  follows  that  3z  =  z 

2  \ 

the  term  ---  3o^  can 

'  z  I 

standard  deviation  is  a 


corie.sponding  to  0  <  z 


j3o^  requires  that  the  change  of  be  much 

of  z  .  Sut  since  z  =  0  at  some  point  near  0 

.  Tha  maximum  z  being  considered  is  Hence 

certainly  be  ignored  if  the  variation  3o_^  of  the 

small  fraction  of  o  over  the  0  interval 
z 

<  o  To  eliminate  the  last  term  of  Eq .  (5a) 

z 


The  argument  has  been  conducted  for  the  finite  interval  0  £  1  o 

To  derive  from  it  the  point-wise  condition  required  for  the  simplification 

of  Eq .  (5a)  one  need  only  postulate  linearity  of  Y(6)  and 

neighborhood  of  0  . 

0 


one  need  only  recognize  that  the  error  function  has  a  maximum  value  of 
unity.  Hence  variation  of  by  only  a  small  percentage  over  the  0  range 
of  interest  makes  the  term  negligible  compared  with  the  term  in  Sz”  for 
z"!  <  o 

It  appears  reasonable — and  numerical  examples  worked  out  confirm 
this  by  on  ample  margin — that  the  variation  of  over  the  operating  6 
range  of  a  functioning  tracker  should  Indeed  be  a  snail  fraction  of  o^. 
Hence  Eq .  (9)  is  a  good  approximation  to  Eq.  (5)  for  all  6  such  that  the 
condition  l^l  ^  is  not  violated  drastically.  Under  these  conditions 
the  6  range  for  which  Eq .  (9)  is  satisfactory  has  a  cumulative  probability 


very  close  to  unity. 


JI.  Conventional  prickers 
A.  General  Relations 

An  elementary  version  of  the  conventional  split  b.an  (phase)  tracker 
is  shown  in  Figure  2.  Suppose  that  the  delay  of  the  signal  from 
hydrophone  to  hydrophone  is  t^  while  the  delay  of  the  interference  from 
hydrophone  to  hydrophone  is  A.  Then 


signal  g^(t)  .  Si[t+(j-l)tJ 

interference  ij(t)  =  lj[t+(j-l)A] 


Figure  2 
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(12) 


If  the  delays  t 


j 


ere  adjusted  so  that 


T 


j 


(j-l)ti 


then  the  suramer  outputs  are 
M  , 


V^(t) 


and 


^|s^[t+(j-l)(t^-tj)]+  ij[t+(j-l)(A-tp]+  nj[t-(j-l)tj|  (13) 

Mr  N 

j(t)  -  Sj|^t+(j+M-l)(t^-tj)j  +  ij[t+(j-rt!-l)(A-tj^  +  nj^j[t-(j-fM-l)tJl 

j-ll  J 


If  tj^  ■  t^  the  array  is  steered  on  target.  The  time  delay  tj^  is  related 
to  steering  angle  6  (6  <>  0  corresponds  to  broadside)  through  the  equation 

t,  «  —  sin  6 
1  c 

d  is  the  spacing  between  hydrophones  and  c  is  the  velocity  of  sound  in  water. 

H^(ja))  and  Hg(ju))  in  Figure  2  are  linear  filters  whose  transfer 
functions  remain  unspecified  for  the  moment.  For  proper  functioning  of 
the  tracker  they  should,  of  course,  differ  in  phase  shift  by  90°  over  the 
entire  frequency  band  processed.  The  symbols  1*^(1)  ^*'*1  will  be  used 

to  designate  the  weighting  functions  (Impulse  responses)  corresponding  to 
the  transfer  functions  H^(ju)  and  Hg(ju)  respectively.  In  terms  of  this 
nomenclature  the  output  y(t)  of  the  multiplier  is 


r 


y(t) 


do 


h  (o) 
A 


v^(t-o) 


►j  dp 


hB(p) 


VB(t-P) 


(15) 


Assuming,  without  loss  of  generality,  that  the  low  pass  filter  has  a 
transfer  function  with  unity  gain  at  zero  frequency,  the  average  tracker 
output  ?  is  given  by 
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Since  x  (t)  and  x_(t)  arc  Gaussian  random  processes,  this  fourfold 
A  D 


average  can  be  expressed  in  terms  of  corrulrtion  functions. 


R  (t)  -  R  (0)  R  (t)  R  (T)  +  R  (t)  R  (-t) 


Recognizing  that 


R  (0)  =  R  (“) 


x.Xb  y 


one  obtains  from  Eq.  (21) 


1  r 


7  I  V *7 

/Ad 

bow  Invoking  Parseval's  theorem 


^  dx  R  (t)  R  (-t) 

^a"-b  Vb 


0  ^  ^  duj  G  (ui)  G  (w)  +  —  /  du) 

j  ""a  "'b  ^  / 


^’x  X 

’‘a  B 


(23) 


(24) 


(25) 


(26) 


Finally,  expressing  the  result  in  terras  of  the  auto-  and  crosscorrclation 
functions  of  v.(t)  and  v„(t), 
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00 

It  J 

1  UU. 

I'^.^dtj)!  [11^(1(11)!  r,  (w)  c  (w)  + 

ll^(.)u))  HB(juj) 

2 

G  (oj) 

V  V 

1 

'  1 

—  03 

'  A  ''b 

^  A  b  ' 

(2 

Equations  (19)  and  (27)  are  the  fundamental  relations  describing  the 
operation  of  the  tracker.  Note  that  they  defend  only  on  the  spectral 
properties  of  the  summer  outputs.  They  will  therefore  hold  equally  well 
when  the  beam-forming  system  contains  provisions  for  nulling  an  interference 
or  for  discriminating  against  an  anisotropic  noise.  The  effect  of  such 


R  with  a  single  subscript  denotes  the  autocorrelation  of  the 
subscripted  variable.  A  double  subscript  indicates  the  crosscorrelation 
between  the  two  subscripted  variables. 


provisions  on  tracker  performance  is  completely  described  by  their 

Influence  on  the  spectral  '■ropeities  of  • 

The  simplest  possible  choice  of  H  (jw)  and  H  (Joj) ,  retaining  only 

A  B 

the  feature  essential  for  tracking , (the  90°  relative  phase  shift)  is 


Hg(ju)  =  1 


(28) 

(29) 


It  will  shortly  become  apparent  that  changes  in  the  filter  functions 
can  be  treated  as  equivalent  changes  in  the  input  spectra  (see  p.  15). 
so  that  Eqs.  (28)  and  (29)  are  now  restrictive  in  any  important  sense. 
Using  Eqs.  (28)  and  (29),  Eq  (19)  becomes 


z  =  2  j  du  ^  (y)J 

0 


V  V 

A  B 


(30)' 


while  Eq.  (27)  reduces  to 

2  2ti 
z  T 


r 

f  1 

2’ 

/  d.j 

G  (m)  G  (m)  -  i 

G  (w)' 

1  V.V-  1 

} 

1 

L  A  ® 

[  A  B  J 

1 

(31) 


Equations  (30)  and  (31)  will  form  the  basis  for  most  computations  in  this 
report . 

B.  Average  Tracker  Output 

According  to  Eq .  (30)  the  average  tracker  output  depends  only  on  the 
cross-spectral  density  of  v  (t)  and  v  (t) .  This  quantity  is  easily 
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computed  from  the  corresponding  autocorrelation  function.  Using  Eqs.  (13) 
and  (lA)  and  assuming  Independence  of  the  ambient  noise  from  hydrophone 
to  hydrophone  (as  well  as  the  independence  of  signal,  Interference,  and 
ambient  noise  from  each  other) 


1 


Oi 


)  stands  for  the  imaginary  part  of  the  bracketed  quantity. 


A-1  1 


R  (t)  =  E 

V  V 

A  B 


H  M 


i=l  k=l  ■ 


X  <Sj^[t+T+(k--M-.\)  (t^-t^)|  +  1 J  t+T+(k+M-l)  (A-t j^)]  +  n^|_^j^['t+T-(k+M-l)tj^ 


M  M 


M  M 


S,=  l  >1  e=l  k=l 

(32)^ 

P,,(t)  and  Pj.(t)  art  the  normalized  autocorrelation  functions  of  signal  and 
interference  respectively.  S  is  the  average  signal  power  and  I  the  average 
interference  power. 

Applying  Eq ,  (17)  to  Eq .  (32)  one  ohta.ins 


M  M 

G  ^  (w)  =  s  )  y  gg(u) 


ju,(k->+M)(t  -L.) 

o  1 


A  B 


i=l  k=l 


M  M 

+  I  ')  8^(w)  e 

d-l  k=l 


jwCk-H+M) (A-tj) 

(33) 


gg(u))  and  gj(w)  are  the  normalized  spectral  densities  of  sigv al  and 
interference  resoectively ,  i.e. 


gg(u.)  Ju  =  j  gj(i.)  du  “  1 

J 

—or* 


It  is  clear  from  Eqs.  (19)  and  (33)  that  modification  of  H^(ja))  and 

lijjtjuj)  by  the  sane  factor  ll(jm)  Is  equivalent  to  modification  of  gg(w) 

2 

•and  g^(j)  by  the  factor  lH(,iu)|  .  As  anticipated  on  p.lA  one  can  therefore 

study  the  effect  of  filtfring  operations,  such  as  prewhitening,  simply  by 

2 

considering  appropriate  nodlf Icat ions  in  the  input  spectra."” 


1  , 
e 


designates  the  expectation  of  the  bracketed  quantity. 


^A  similar  statement  holds  concerning  the  effect  of  filters  on  the 
output  fluctuation  [see  Eq .  (31)]. 
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I 

I 


Substitution  of  Eq.  (J3)  into  Eq.  (30)  yields 

CO 

M  M 


z  =  2^  \  I  du  Ssf!  (u>)  sinfuiCk-i+M)  (t  -t  )j  +  1  f’.Cui)  sin fioCk-il+M)  (A-t  )  j 

L-i  C-i  I  S  ol-^  I  1 

11=1  k=l  1- 

(34) 

Equation  (34)  will  now  be  evaluated  for  var  .ous  signal  and  interference 
spectra. 

1)  Uliite  Spectra.  Consider  first  the  case  of  sign^^J  and  interference 
spectra  whitf  (or  prcuhitened)  over  the  entire  processed  frequency  band. 

Then 


‘  I 


8s(‘")  =  =1 


1 

0 


UJ  Id 

I  I  /. 


'  u  >  u 

■  »  -  Q 


(35) 


Equation  (34)  is  now  easily  integrated,  with  the  result 

< — I  *  ”  cos  [(k-ili-M)^  (t  -t  )J  1  -  cos|(k-{+M)  w  (A-t,)]  I 

7=  )  \  Js - ^ - ^ - 2 - 

oTi  Jl",  (k-?+M)  0)  (t  -t.)  (k-i+M)  u  (A-t.) 

i- 1  K-  1  I  O  O  1  O  1 


Since  the  indices  k  and  H  occur  only  in  the  conbination  k-Jl,  it  is 
cor  .ent  to  introduce  the  notation 

k  -  il  =  r  (37) 

Then  Eq .  (36)  becomes 


M-1 


1 


cos  [(r+M)u^(t^-t^)  J  1  -  cos[(r+M){ii^(A-tj)] 


+  I 


r=-(M-l)| 


Introducing  -he  notation 


(r+M)u)  (A-t.) 
o  1 


x  =  Id  (t  -t  )  and  y,  =  w  (A-t  ) 

1  o  o  1  1  o  V 


(M-  lr|) 


(38) 


(38a) 


one  can  rewrite  Eq .  (38)  in  Che  forn 


z  -  S 


1  -  cos(r+M)x,  T  1  -  cosRr-m)  (x,+y,)l 

\  < - L+  I  - L - i — LlL^fM-Irl) 

Lj  (r-m)x,  S  (r+M)(x,4y  ) 

_  /\J  1  \  A  A  ^ 


r— (M-1) 


(39) 


Here  x^  is  a  measure  of  the  array  steering  angle  relative  to  the  target 
bearing  while  is  a  measure  of  the  Interference  bearing  relative  to  the 
target  bearing.  The  relation  of  these  two  quantities  to  the  actual  target 
bearing  9^,  Interference  bearing  9^,  and  steering  angle  9  Ig  specified  by 


the  expressions 


t  ■  sin  9 
o  c  o 


A  «  —  sin  9t 
c  I 

t,  ■  —  sin  9 
1  c 


(40) 

(41) 

(42) 


Figure  3  shows  Eq.  (39),  normalized  with  respect  to  S,  plotted  as  a 
function  of  Xj^  for  j  ■  5,  M  =  20  (40  element  array),  and  y^  «  -5.^  With 
a  broadside  target,  *  2ir  x  5000  rad /sec,  and  1  ft  hydrophone  spacing, 
yj  ■  -5  corresponds  to  an  Interference  bearing  of  about  53°  from  the 
target.  Figure  4  gives  a  similar  curve  for  yj^  ■*  -2,  corresponding  to  an 
interference  bearing  of  about  19°  under  the  same  assumptions.  Figures  5 
and  6  give  equivalent  curves  for  M  =  10  (20  element  array) .  In  each  case 
the  plot  exhibits  generally  the  expected  form  (see  Figure  1)  near  the 
target  bearing.  It  shows  a  similar  functional  behavior,  greater  In 
amplitude,  near  the  Interference  bearing.  Perhaps  the  most  striking 
feature  of  the  curves  is  the  offset  in  the  axis  crossing  caused  by  an 
interference  only  7  db  above  the  signal  level,  even  when  the  interference 
is  relatively  remote  in  angle  from  the  target.  If  the  tracker  indicates 
a  bearing  corresponding  to  the  zero  of  Y,  then  this  effect  can  lead  to 
appreciable  systematic  error  in  indicated  bearing.  With  M  ■=  10,  y  -2 


^The  curve  on  the  right  gl res  the  complete  pattern,  while  that  on  the 
left  shows  the  neighborhood  of  che  origin  In  expanded  form  (target  response). 
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the  zero  occurs  at  «  0.078  .  Vlth  a  broadside  target,  >  2ir  x  5000 

rad/sec,  and  d  ■  1  ft,  this  is  equivalent  to  c  systeoatlc  error  of  0.71*^. 

Even  with  H  •  20  and  y  >  -5  the  systematic  error  is  still  about  0.06°. 

The  importance  of  the  problem  clearly  depends  on  the  relative  strength 

of  target  and  interference.  To  exhibit  this  effect  more  clearly.  Figure  7 

shows  the  offset  in  the  zero  (systematic  error)  as  a  function  of  for 

M  «  10  and  M  ■  20  and  various  values  of  y^ .  With  the  numerical  values 

of  and  d  used  previously,  the  vertical  scale  can  be  converted  to  degrees 

57  3 

of  systematic  error  through  multiplication  by  - j — ' - *9.1  .  The 

0)  —  cos  e 
o  c  o 

curves  for  y^  ■  -I  and  y^  «  -2  terminate  within  the  range  of  the  graph 
because  larger  values  of  -I-  cause  the  zero  near  the  true  target  bearing  to 
disappear  entirely.  An  instrumentation  attempting  to  Crack  this  zero 
would  therefore  fall. 

2)  Algebraic  Spectra.  The  relatively  large  systematic  errors 
introduced  in  the  above  calculations  by  rather  remote  sources  of  Interference 
are  due  at  least  in  part  to  the  slow  decay  of  Eq.  (39)  away  from  x  ■  0  . 

That  effect,  in  turn,  may  be  attributed  to  Che  sharp  decay  of  gg(u)  and 
g^(u)  at  One  suspects  that  gradually  decaying  spectral  functions  would 

produce  more  favorable  results.  More  precisely,  since  the  Integrals  in 
Eq.  (34)  are  closely  related  to  Fourier  transforms  of  gg(u),  less 
systematic  error  can  be  expected  from  spectra  whose  Fourier  transforms 
(autocorrelation  functions)  decay  rapidly  away  from  the  origin.  An 
analytically  convenient  function  with  these  general  properties  Is 

gg(u)  -  gj(w)  -  - 2  (43)^ 

(u)^  +  U)  ^) 

a 

^Note  that  this  spectral  function  falls  to  zero  at  u  ■  0  ,  a  feature 
generally  present  in  practical  systems. 
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Substituting  Eq.  (43)  into  Eq.  (34)  one  obtains 


Z  m  2U 


^  /  do) - - 2  sin  [u(k-t.+M)  (t^-t  j)]  +  I  sln[<ji)(k-l+M)  (A-t  j)j| 

1^1 -i  (u^+u;  h  ^  J 


M  Mr 


■  fE  E^ 


.S(k-«,+M)a)  (t  -t.)  e 
I  a  o  1 

1-1  k»l^ 


-|(k-14M)(t^-tj)a)^| 


+  I(k-l-«l)u  (A-t,)  e 
a  1 


Now  changing  the  index,  as  before, 

k  -  t  «  r 


-|(k-)l4M)a)  (A-t.) 

a  1 


M-1  I' 

i  E  <S(r+M)a,^(t^-t^)  e 
r=-(M-l) L 


-1 (r4M)(t  -t  )w  1 
Ola 

-l(r+M)(A-t  )a)  n 

+  I(r+M)(A-t,)a)  e  ®  MM  -  r  ) 

la  I 


(44) 


(45) 


(46) 


Finally,  Introducing  the  notation 


a  o  1 


u  (A-t  ) 
a  o 


(47) 


one  obtains 


z  ■ 


It  r  -|(r+M)x  I 

Y  ^  ^S(r+M)x2  e 

r— (M-1)  I 


+  I(r+M)(x2+y2)  e 


-| (r+M) (x2+y2)l 


>(M-  |r|) 


(48) 


Figures  8  and  9  present  plots  of  z*  equivalent  to  Figures  3  and  6.  In  both 
cases  the  plots  show  a  much  smaller  Influence  of  the  Interference  on  the 
behavior  of  ?  near  the  target  bearing  than  was  the  case  with  white  spectra. 
Some  caution  must  be  used  in  comparing  the  systematic  errors  directly, 


A-21 


Form  EC'12 


DUNHAM  LABORATORY,  YALE  UNIVERSITY 


A-22 


because  the  physical  meaning  of  is  not  exactly  the  same  as  that  of 

Since  a  substantial  amount  of  the  total  power  of  the  algebraic  spectrum 

lies  above  u  ,  a  fair  comparison  clearly  demands  that  u  be  chosen  smaller 

than  u^.  One  basis  of  comparison  with  a  certain  intuitive  appeal  is  the 

following;  Suppose  that  the  spectrum  is  shaped  primarily  by  the  filters 

H.(Ju)  and  H  (ju).  Then  it  is  reasonable  to  assume  that  the  ambient  noise 

has  the  same  spectral  form  as  the  signal.  We  choose  w  such  that  the 

A 

detection  index^  for  both  types  of  spectra  is  the  same  in  the  absence  of 
Interference.  A  straightforward  computation  yields  * 

Figure  10  presents  information  analogous  to  Figure  7,  the  systematic 
error  for  the  case  of  algebralr  spectra,  plotted  as  a  function  of  for 

g 

various  values  of  y,  *  ~  y^  •  Also  shown  for  the  sake  of  comparison  are 

i  Tf  Z 

two  of  the  curves  from  Figure  7  (white  spectra) .  The  vertical  scale  is 

g 

the  displacement  of  the  zero  measured  in  units  of  x^  ■  —  X2  .  Hence  for 

u  >  2n  X  5000  rad/sec  (w  »  2n  x  I960  rad/scc)  and  d  «  1  ft  ,  the  vertical 
o  a 

scale  must  be  multiplied  by  9.1  to  obtain  systematic  error  in  degrees 
(exactly  as  with  Figure  7) . 

3)  Exponential  Spectra.  Another  analytically  convenient  spectral 
function,  intermediate  in  its  cutoff  properties  between  those  of  Eqs.  (35) 
and  (43) ,  is 

gg(u)  ■  gj(w)  =  e  ^  (49) 

b 

The  resulting  expression  for  T  is 


z  ■ 


M^l  z' 


2  2 
(r+M)^X3‘' 


+  I 


(Xl+Fl) 


1  +  (r+M)^(x3+y3)^  j 


f  (M-|r|)(M  +  r) 


(50) 


^Report  No.  3,  Eq.  (49). 
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where 

*3  “  ^3  " 

Figure  11  shows  the  systematic  error  plot  corresponding  to  Figures  7 

and  10.  For  proper  comparison  under  a  criterion  requiring  equal  detection 

(1) 

Index  In  the  absence  of  Interference,  one  finds  here  that 

expected,  the  systematic  error  Is  Intermediate  to  that  produced  by  the 

other  two  spectra.^ 

C.  Output  Fluctuations 

According  to  Eq.  (31)  one  requires  the  spectra  G  (u)  and  G  (u)  In 
addition  to  the  already  computed  G  (u)  [Eq.  (33) j  .  Computing  first 
the  autocorrelation  functions,  one  obtains  from  Eqs.  (13)  and  (14) 


A  B 

MM  MM 

-  ^  Pg[T+(k-l)(t^-tp]  +  +  M  N  0k(t) 

k-1  1-1  k-1  1-1 

(52) 


P^(t)  Is  the  normalized  autocorrelation  of  the  ambient  noise  and  N  Is  the 
average  ambient  noise  power. 

Fourier  transforming  to  obtain  a  corresponding  expression  In  terms  of 
spectral  functions. 


M  M 

G  (w)  -  G  (u)  -  sV  )  g_(u)  e 
*  ®  k-1  1-1 

M  M 


J(k-l)(t^-tj)u) 


^  j(k-l)(A-t,)(.i 

+  12^  gj(<*))  e  ^  +  M  N  gjj 


(w) 


k-1  1-1 

Inserting  Eqs.  (33)  and  (53)  Into  Eq.  (31)  and  using  the  index 

r  =  k  -  1 


(53) 


(54) 


Hiote,  however, 
terminate  at  smaller 


that  the  curves  for  exponential  spectre  consistently 
I/S  than  those  for  the  other  spectral  types. 
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one  obtains 


2  2?i 


1^  r=-(M-l) 


.ir(t^-ti)c. 


(M-lr|) 


M-1  ,  /  N 

r-  Jr(A-t  )a. 

+  I  )  Pj^(w)  t-  (M-lr|)  +  M  N  gj^Cei) 

r='-(K-i  ) 


^  jiw(r+M)(t  -t  ) 

/  ggCiiOe  “  (N-|r|)+I  ^  gj(w)e 


jui(r+M)  (A-t  ) 


(M- 1 r I ) 


r=-(!;-L) 


r=-(M-l) 


In  tile  cases  of  prlninry  interest  here  the  signal  power  is  much  smaller 
rh.an  the'  int  rfirenci  power,  (S  may  or  may  not  be  small  conpared  with  N.) 
In  tint  ease  the  slgnil  contribution  to  the  output  vari.ni'Je  is  negligible 
a^d  ifi .  (ii)  reduces  to 


■  2  _  2r  2 

z  “  T  ^ 


M-1  M-1 

r=-(M-  1)  ci  =  -(M-l) 


).e(/+(|l  (A  t  )  jiiifiM  q  I  2MU.'’i-t  , ) 


(M-|r|)  (M-|q|) 


l'r(A-t^)  a  2  2  I 

+  ^  ‘  ^  gj(u)  ,H-|r|)  +  M“N  gj^,  (u) 


r--  (ii- 1 ) 


Lquation  (56)  must  now  be  evaluated  for  the  various  spectra  considered 


previously . 


1)  White  Spectra.  Consider  first  the  speetrai  functions 


gj(w)  “ 


I  c  I  2  tv., 


Substitution  into  F.q.  (56)  yields 


i'°  I  ,  J(,.>(r+q)(A-t  )  +Jtj''r+q  +  2H)(u-t  ) 

|  d.  e  ^  -c  '  J(M-|ri)(M-|ql) 

^*^0  r*=-(M-l)  q“-(M-l) 


+  2  M  N  1 


j.r(A-t.) 


^  (M-  I  r  I )  +  M 


Ml  Ml  sin|(r+q)(A-t  )ui  ]  sin  [( r+q+2M)  ( A- 1  )  u  j 

»  r.-(M-l)  1  ° 


+  2  M  N  I 


r=-(M-l) 


sin  I  r  (A-t,  )ui  1 
1..0. 

r  (  A  - 1 ,  )  w 
1  o 


The  symbols  x  and  y  are  defined  in  Eq.  (38a).  x  +y  =  (A-t  )uj  is 
i  1  1  1  1  0 

clearly  a  measure  of  steering  anglu  relative  to  the  interference  bearing. 

Figure  12  shows  a  normalized  version  of  o^,  the  square  root  of  Eq .  (59), 

N 

plotted  as  a  function  of  Xj^+y^^for  M  =  10.,  20  and  several  values  of  y  •  De¬ 
pendence  on  the  latter  parameter  is  small,  as  one  would  expect  in  an 
interference- 


dominated  environment.  If  the  tracker  Is  to  function  properly,  o  must 

z 

certainly  be  small  compared  with  the  peak  value  of  z  near  the  target 
bearing.  According  to  Figures  3  and  5  Tiorizontal  excursions  must 
therefore  be  confined  to  a  range  small  compared  with  0.1  and  0.2  for 
M  ■  20  and  10  respectively.  Over  such  a  range  Figure  12  exhibits  changes 
equal  to  no  more  than  a  small  fraction  of  the  value  of  o  at  any  point 
within  the  range.  Thus  the  approximation  leading  to  Eq.  (9)  Is  amply 
justified.^ 

For  remote  Interference  Eq.  (59)  reduces  to 


2  _£ 
3  T<i) 


{■ 


I  <M(2M  +1) 


+  3  Y 


(60) 


If  M  »  1  one  can  therefore  write  approximately 


3  T(d 


3  2 

M-*  I 


for 


f- 


(61) 


Much  as  in  the  case  of  detection  one  finds  the  fluctuation  power  varying 

3  2 

with  M  In  an  interference-dominated  environment  and  with  M  In  an  amblent- 

nolse-dominated  environment.  The  dividing  line  between  the  two,  again 

as  In  the  case  of  detection,  depends  on  s'  t  even  a  relatively 

weak  Interference  can  dominate  tlie  output  fluctuation  if  the  array  Is 

sufficiently  large. 

2)  Algebraic  Spectra.  Next  consider  again  the  algebraic  spectra 

gj(u)  ■  - 2 

(u)^  +  m^) 

^The  approximation  Is  poorest  for  very  small  This  is  the 

condition  of  steering  angle  very  close  to  Interference  bearing,  where  no 
satisfactory  tracking  performance  would  be  expected  In  any  case. 
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Substitution  Into  hq .  (56)  yields,  after  the  usual  algebraic 
simplif icatlons  , 


2  2  -|(!'+q-2M)(x,+y,)  I  ,  . 

z  sTiT' ^  ^  1+1  (r+q-2M)  (x2+y2)  i-  -jl  (r+q-2M)  (x^+y^)'! 


Z_,  ^ 

[r=l  q=l  !| 

-I  (r+q) (x^+y^)  I  r 


J  +  1 (r+q) (x  +y^) I  -  ( r+q ) (x  +y  ) | ^  I  (M- | r-t1  ] ) (M- | q-H | ) 


i-  2M  j  11+2  ^  e 


y-  -rlx^+y,!  I 


l+lr(x2+y2)|  -  "3 1  ^  ^  (M-r)  +  “2 


For  1X2+22!  ^  (steering  angle  remote  from  interference)  and 

-I  =  77  u.  t'.q .  (63)  becomes  identical  with  Eq .  (60). 
a  o  o 

N 

FUot.'-.  of  Eq.  (63)  for  several  values  of  M  and  j  are  shown  in  Figure  13, 
Fere-,  as  in  Figure  12,  the  approximation  leading  to  Eq .  (9)  if  c’early 
satisfied . 

D.  Random  bearing  Error 

According  to  Fq.  (9)  the  random  component  of  bearing  error  in 


evaluated  at  the  target  bearing.  Only  the  slope  of  the  average  respon.se 
curve  (tlie  sensitivity)  remains  to  be  calculatc^d .  For  w'nite  spectra  one 
obtains  from  L,q .  (39) 

3x 

j_z  _  3  z 1 

v'  “  3Xj^'  30 


Su.  —  cos  0  \  -- 
o  c  c\  2 


!.sln(M+r)y  l-cos(M+r)y  ( 

2-«  *5  ^  |-(TiT-07---  ,2  ^ 

r=-(M-l)l-  I 


A-;!2 


l  orm  r.E*l2 
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A  similar  computation  for  the  case  of  algebraic  spectra  yields 


3  z 
30 


3x 

3z  2 
3x^  30 


M-1 


rSu 


cos  u 


-  M^+  — 
2  2S 


l-(M+r) ty„ 


-(M+r) |y 


2  I  ,  .  \ 

(t!-lr  |)  (M+r)  / 


(M-l) 


(06) 


For  remoti!  interference  and  u  =  tt  •  bps.  (65)  and  (66)  are 

a  0  0' 

related  by 

3z  I  _  37, 1 

30  I  “  '  8  30’ 

ialgeb.  white 


(67) 


Normalized  plots  of  Eqs.  (65)  and  (66)  for  various  values  of  M  and  — 


are  shown  in  Figures  lO  and  15  respectively.  In  each  case,  bvit  particularly 
fur  tile  algebraic  spectrum,  the  sensitivity  quickly  .ipproaches  its 
asymptotic  value  once  the  separation  of  interf ertnet-  and  target  exficds 
the  basic  beamwldth  of  the  array. 

Normalized  plots  of  random  bearing  error  [Eq.  (66) J  are  shown  i:. 

N 

Figuies  16  and  17.  S.nci  the  -'aiucs  of  y  +'■1  correspond  to  an  interference- 
dominated  envlronirent ,  the  dependence  on  this  parameter  is  quite  weak.^ 

Fince  S/I  appear.s  oi. '  y  in  the  expression  for  sensitivity  it  is  clear  from 
Fig, ores  lA  and  15  that  this  parameter  has  little-  influence  on  the  random 
bearing  er-or  except  for  inter fere-nces  within  .a  bi  am  width  of  the  target  - 
in  whicli  cas'..  -liscussion  of  the  fluctu.ation  error  becomes  an  academic  matter 
because  the  postc.lated  tracker  can  no  longer  distinguish  hi.t\,'een  target  and 
interference.  Conversion  of  the-  vertical  scale  of  Figure  16  into  degree's 
of  rms  error  1,  .accomplishi  d  through  inu  1 1 ipl icat i  r.  hy 


57.3 


I 

S 


_J _ 

d 

—  cos  '■ 


(68) 


Hiote  =  2  ,  r  =  10  is  fairly  close  to  tli.  borderline  of 

imbiv.nt  noisi  domination,  hence  the  deviation  i,s  great-  si  in  this  casi.  . 
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Of  considerable  pmctic::!  inportance  is  the  question  whether  the  random 
component  of  error  cnrrvntly  under  di:;cussion  is  larger  or  smaller  than  the 
systematic  error  investigated  in  Section  II. B.  In  order  to  avoid  tedious 
comparison  of  a  1  -.rge  nut:ber  of  .ta.scs,  the  following  a.pproximatc  analysis 
is  carried  out: 

If  tne  tracker  is  to  function  properly  the  shift  in  the  null  (systcnatlc 
error)  should  be  confined  to  a  region  sufficiently  close  to  the  true  bearing 
so  that  th  average  bearing  response  pattern  (Figure  1)  may  be  regarded  as 
linear  Dver  the  region  in  question.  In  that  case  one  can  write 


Svster  'itjc  ertor  (in  r'uiians)  = 


(69)' 


3  z 

ae 


io 


where  z^  is  the  on-target  value  of  tin  average  b'.-arlng  response.  Hence 


from  Eqs.  (6A)  and  (69) 


z 

systematic  error  _  _ o 

rms  r.andoc  trror  z 

o 


(70) 


Is  obt.iined  by  tting  =  0  in  F.q .  (39)  and  x^  =  0  in  Eq .  (A'),  . 

r.  iiTill  arly  "quatlnns  (S'')  .nd  (<'3)  yield  ■■  .  Figures  )  3  and  19  show 

o 

normalized  plot.s  of  Eq .  (70)  on  a  logarithmic  scale  ior  whltt'  .'.nd  algebraic 
spectra  respectively.  An  Inter  fen nco-doninated  envlronr’ent  is  assumed, 
so  that  dependence  of  the  noise  level  on  N  may  be  igno'ed-  F.xcept  for 
the  neighborhood  of  y^  =  ?.t  tlu^  curves  for  whltt  sp  ctra  (Figure  18) 
are  essentially  straight  linos  with  a  slope  of  ~2^  >  once  the  interference 
is  separated  from,  the  target  by  mor'  tlian  the  width  of  the  bear,  pattern." 
This  behcavlor  has  important  practlc  -.l  impllc.ations .  The  abscissa  y^^  is 
related  to  target  .and  inti  rferenco  bi  .aring  through  the  equation 


Ihij  is  clearly  equivalent  o  the  .-ipproxim.’.tlon  of  Equation  (9). 

Each  of  the  curves  exhibits  a  sharp  downward  excursion  near  =  2n  . 
To  avoid  confusion  on  ehe  graph  this  has  been  show'n  only  for  M  =  5  . 
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(71) 


If  (u^'  =  Ku)^  ,  the  corresponding  is  Ky^^.  Because  of  the  (-  y)  slope 

-1/2 

this  results  In  R'  =  K  R  .  Hence 


so  that  the  error  ratio  remains  unchanged.  In  other  words,  the  error  ratio 
is  independent  of  the  cutoff  frequency  as  long  as  y^  lies  on  the  linear 

portion  of  the  curves.  Furthermore,  note  that  the  horizontal  spacing  between 
curves  is  essentially  constant  at  a  value  corresponding  to  a  factor  of  2. 

An  increase  in  the  number  of  hydrophones  decreases  the  error  ratio  if  all 
other  parameters  remain  fixed.  This  implies  that  the  total  array  length 
increases  linearly  with  M.  In  practice,  it  may  be  more  realistic  to 
consider  the  array  length  fixed  and  to  vary  the  hydrophone  lacing  d 
inversely  with  M.  Assuming  that  phone-to-phone  independence  of  the  ambient 
noise  is  maintained  during  these  changes  one  finds  that  to  double  M  implies 
cutting  y^^  in  half.  Hence  the  error  ratio  is  unaffected  by  the  change. 

One  concludes  that  for  fixed  overall  array  length  the  i-rror  ratio  is 
independent  of  the  number  of  hydrophones  as  long  as  y^^  lies  on  the  linear 

portion  of  the  curves . 

The  next  question  is  clearly  whether  y^  does  or  does  not  lie  on  the 
linear  portion  of  the  curves  in  situations  of  practical  interest.  One- 
need  rot  be  greatly  concerned  with  the  nonlinearity  for  small  y^,  for  one 


A-  12 


would  hardly  expuct  to  operate  the  system  with  interferences  separated  from 

the  target  by  less  than  a  beam  width.  To  obtain  an  upper  bound  on  y^ 

consider  d  ••  1,  u  >  2v  x  5000  rad/sec.  With  a  broadside  target  (6_  *  0) 
o  T 

£q.  (71)  yields  a  maximum  y^  of  2ir  .  Thus  all  but  a  small  range  of 
interference  bearings  near  endfire  produces  values  of  y^  within  the  linear 
portion  of  Figure  18. 

It  is  now  possible  to  determine  the  value  of  the  error  ratio  in  the 
linear  region,  /hich  is  independent  of  all  parameters  except  overall 
array  length,  smoothing  time,  and  bearing  angles. 


systematic  error 
rms  random  error 


.  53  Vf- 


sin  6-  -  sin  6 
I  o' 


where  L  is  Che  array  length  in  feet.  For  smoothing  times  of  Che  order  of 
seconds  and  moderate  array  length  this  ratio  is  clearly  larger  than  unity. 

The  error  ratio  curves  for  algebraic  spectra  [Figure  (19)]  do  not 
exhibit  all  of  the  invariance  properties  present  in  the  white  noise  case. 
Confining  attention  to  the  region  y2  S  ^  (corresponding  to  y^  >  '^)  one 
still  finds  a  horizontal  spacing  roughly  proportional  to  M  and  hence  approxi¬ 
mate  invariance  of  the  error  ratio  relative  to  M  (for  fixed  array  length) . 
However,  the  slope  of  Che  curves  below  *  2  exceeds  ^  .  The  reasoning 
of  Equations  (71)  -  (73)  therefore  leads  to  Che  conclusion  that  the  error 
ratio  is  a  monotone  decreasing  function  of  .  The  conversion  from  Che 
normalized  R  plotted  in  Figure  19  Co  actual  error  ratio  is  accomplished 
through  Che  relation: 


systematic  error  _  2/T 
rms  random  error  ir  a 


(75) 


With  a  40  ft.  array,  1  ft.  hydrophone  spacing,  and  u  “  10,000  rad/sec 

3 
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this  leads  to 


^t-ema^.ic  error  ^  ^  ^  ^  2 

rms  random  error  —  - 


(76) 


Thus  the  error  ratio  is  st^ll  in  excess  of  unity  for  most  reasonable 
smoothing  times.  Only  with  larger  arrays  or  wider  bandwidth  is  the  random 
error  likely  to  become,  dominant. 


I  1 


rep'^escriLs  split  version  of  the  simples;  conf ipuration  analysed  in 

Report  No.  /II  (■■•.  e  ll,  ure  1  ot  that  report).  The  delays  are  adjusted 

In  accordar.fc  with  tlu  Intt- :  f  erence  delay  from  hydrophone  to  hydrophone, 


so  that  the  subsequent  subtractions  achieve  In  principle  complete  elimina¬ 


tion  of  the  Interference.  The  delays  6^  align  the  remaining  signal  com¬ 
ponents.  From  this  point  on  the  Instrumentation  is  Identical  with  that  of 
the  conventional  tracker  (Figure  2).  Thus,  as  pointed  out  previously, 
Fquations  (39)  and  (27)  |  or,  with  the  previous  assumptions  concerning  H 

and  Hg  ,  Equations  (30)  and  (31) J  remain  valid.  The  spectral  functions 
G  (ui)  ,  G  (oj)  aiid  G  (qi)  must  now  be  computed. 

From  Report  No.  21,  Equation  (27)  : 

M-1  ... 

K“i 


and  V  (t)  “  y  \  8  J  t  +  (k-Jffi)  (t  -t,)|  -  a  J  t  +  (Ic+m)  (t  +^)  i 

^  OXJ  lu  OxO  '< 


+  Vk-f-l/  -  ^  ^ 


It  follows  that 


M-1  M-1 


R  tT)  -  R 
A 


-  p  |T+(k-i)  (t  -t.)  -(t  -A) 
H  L  o  1  o 


+  n[2(M-1)  Pj^(t)  -  (H-2)  [t+Ctj^-A)]  -(M-2)  P^|_T-(t^-A)  ]  j 


Fourier  transforming; 


G  (u)  ■=  G  (m)  “  S  g  (aO  )  j?-  (to-tj^)  +  (to-A)J 

''a  ''b  ®  k^^l  jS'l  ^ 


gjojUk-1)  Cto-til-Uo-A)] 


,1w(t  -A) 


/  J  W\  ,  «J/ 

+  N  gj^(o)  |2(:i-l)  -(■•.-2)  e  -(M-2) 


-j(j(t  -A) 


Furthermore , 


M-1  M-1 Jr  1 

(t)  -  S  V  )  (2p  T+(M+k-l)  (t  -t.) 

B  °  ^-1 

Ps|^T+(M+k-4)  (tp-tj)  +  (tp-a)J  -  |^T+(M+k-t)  (tp-tj^)  - 


1 


(81) 


Hence 


G  (o))  “  S  g  (dj) 
v»v_  , 

A  B  k==l  t-l 


\ 

_  gja)i(IH-k-0  (Vtj^)  -  (VA)j|  (82) 


As  in  Che  case  of  the  conventional  tracker,  it  is  clear  from  Equations 
(19),  (27),  (80)  and  (82)  that  multiplication  of  h^(ju)  and  Kg(j(o)  by 
K(ju))  Is  entirely  equivalent  to  multiplication  of  gg(w) 

H(Jii))|^  . 

Substituting  Equation  (82)  Into  Equation  (30)  one  obtains  Che  ex¬ 
pression  for  the  average  tracker  output 


_  M-1  M=,l  ,a>  j 

*  -  2S  r  do,  g  (w)  ( 

k-ll-lio  I 


2  sin  (i)(M+k-l)  (t  -t,) 
o  X 


-  sin  uj^(M4-k-l)  (t^-tj^)  +  (tg-A)  -  sin  ioj^(M+k-ll)  (tj,~tj^)  -  (t^-A)j  | 


with  Che  change  of  index  k-1  ■  r  this  becomes 


z  -  2S 


M-2  «  /  r  1 

)  du  g„(ii))  (2  sin  (i)(r-fM)  (t  -t,) 

r-TM-2)>/o  ®  I  L  o  IJ 


sin  (ijKr+M)  (t  -t.)  +  (t  -A)  -  sin  u  (r+M)  (t  -t,)  -  (t  -A)  i  (M-l-|r|) 

L  oj  L  °J| 


(84) 


The  corresponding  expression  for  the  output  variance  Is  obtained  formally 
by  substituting  Equations  (80)  and  (82)  Into  Equation  (31).  However,  in 
the  cases  of  primary  Interest  here  the  signal  power  will  be  small  compared 
with  the  ambient  noise  so  that  one  can  Ignore  the  signal  contribution  to 
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Thfcr  the  variance  of  z  becomes 


the  output  fluctuation.' 


2  2tiN‘'  2  ’ 

0  ^  /  dm  p,.  (m)  ;  2(M-1)  -  (M-2)  e 

^  ^  (  I  2 

^  -jm(t.-A)r 

-  (M-2)  L  ^  ' 


Equations  (8A)  and  (85)  must  now  be  evaluated  for  the  various  spectral 
functions  under  study. 

B .  Average  Tracker  Output 

Consider  first  the  white  spectrum  specified  by  Equation  (35). 


Integration  of  Equation  (8A)  yields  after  some  algebraic  simplification 

11-2  I  1-cos (M+r)  X  1-cos f (M+r)  x  -y  1 
=  S  -2 - ^  - - - ± — — 


r=-TM-2)  \  \ 


(M+r)  x^-y^ 


1-cos  [(M+r)  ^ 

(M+r)  Xj+yj^  I 


(M-l-Iri)  (8ft) 


Xj^  and  y^  arc  defined  in  Equation  (38a). 

The  equivalent  expression  for  the  algebraic  spectrum  I  Equation  (43) 


-  '•W  j  -|(M+r)  x.l  -I  (M+r)  x.-Vjl 

Z  =  Y  ^  2(M+r)  X-  e  -  j  (M+r)  x^-V-J  <=■ 


r=-(M-2) 


-  (M+r)  X  +y  e 


■|(i;+r)  x,,+y2' 


(M-l-|r  I) 


x^  and  y^  are  defined  in  Equation  (47). 

Plots  of  Equations  (86)  and  (87)  for  selected  parameter  values  arc.: 
given  in  Figures  (21''  and  (22)  respectively.  Only  tlu  neighborhood  of 
the  target  is  shown,  because  the  average  output  elsewhere  is  too  small 
to  be  significant  on  the  scales  of  the  two  grapl.s.  The  most  important 


With  S  >  N  thi're  is  presumably  very  little  difficulty  In  tracking 
once  the  interfarence  has  been  elimlr.iited . 
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aspect  of  the  curves  Is  the  complete  absence  of  any  offset  In  the  null 
(systematic  error) ,  a  conclusion  also  apparent  from  the  equations  upon 
setting  -  0  In  (86)  and  =0  In  (87).  In  comparing  Figures  (21) 

and  (22)  one  must,  of  course,  keep  in  mind  that  >’2  ~  (  equal 

detection  Index,  so  that  y2  ■  -  2  corresponds  very  nearly  to  y^^  »  -  5  . 


C.  Output  Fluctuations 

Substitution  of  the  white  noise  spectrum  |Equatlon  (57)j  into 
Equation  (85)  yields  after  algebraic  simplification 


AitK 

Tu 


f 


(M-1)^  +  !s(M-2)^ 


2(M-1)  (M-2) 


sin(Xj^+yj^) 


_  sin  2(x  4y  ) 
+  S5(M-2)  --  -  ^ 


2(Xj^+yj^) 


The  equivalent  expression  for  the  algebraic  spectrum  is 


(88) 


2  2  -|x2‘^2l 

1)  +  h(.K-2)  -  2(M-1)  (M-2)  e 


[1+1X2+72!  - 


3I  2  ^  X2'fy 2 1  r  ,  . 

'’I  +  i5(M-2)^  e  [l+2|x2+y2l 


-  31^2+72 


^1} 


(89) 


Figures  (23)  and  (24)  show  the  standard  deviation  plotted  in 

normalized  form  as  a  function  of  (Xj^+y^^)  and  (x2+y2)  respectively. 
(Xi+yi)  and  (x2+y2)  are,  of  course,  measures  of  steering  angle  relative 
to  the  interference  bearing. 


D.  Random  Bearing  Error 

The  sensitivity  of  the  tracker  is  obtained  by  differentiating 
Equations  (86)  and  (87)  with  respect  to  9  and  evaluating  the  deriva¬ 
tives  at  the  target  bearing.  The  result  is: 
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For  the  white  spectrum 


U)  S  —  cos  0  1 1-2 - +  2  - r - 1 

O  C  ot  y  V  ^ 


M(M-1)‘ 


(90) 


For  the  algebraic  spectrum 

d  r  .  .  "I^ah  2 

-  nu)  S  -  cos  0  1-(1-  y,  )  e  M(M-l)  (91) 

a  c  o  L  2' 

o 

For  remote  interference  and  u  »  (■t/8)  u  (the  basis  of  comparison 

s  o 

used  before)  Equations  (90)  and  (91)  are  related  by 


algeb. 


2 

ii_ 

8 


white 


(92) 


exactly  as  in  the  case  of  the  conventional  tracker.  For  bearings  renote 
from  the  interference  and  oj  =  (71/8)0)  ,  a  approaches  the  same  value 

B.  0  2 

for  both  types  of  spectra  |^see  Equations  (88)  and  (89)  .  Hence  the 
asymptotic  ratio  of  fluctuation  errors  is 


rms  error  with  algebraic  spectrum 
rms  error  with  white  spectrum 


"2  for  remote  interference 

TT 


Plots  of  normalized  rms  error  versus  interference  bearing  relative  to 
target  bearing  (measured  by  y^^  and  y2)  are  given  in  Figures  25  and 
26.  The  qualitative  differences  between  the  two  sets  of  curves  is  due 
in  large  measure  to  the  very  rapid  approach  to  the  asymptotic  value  of 
the  sensitivity  for  algebraic  spectra.  In  other  words,  once  the  inter¬ 
ference  is  separated  from  the  target  by  more  than  a  small  minimum  angle, 
the  sensitivity  at  the  target  bearing  is  almost  independent  of  the  in¬ 
terference  bearing  (for  algebraic  spectra).  This  is,  of  course,  precisely 
the  lack  of  influence  of  the  interference  on  average  response  which 
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motivated  the  Investl station  of  algebraic  spectra  in  the  first  place. 


As  a  final  point  of  interest,  it  nay  be  u.seful  to  compare  the  rms 
error  of  the  conventional  tracker  vltii  tliat  of  the  nulling  tracker. 

From  Equations  (59),  (65),  (88)  and  (90)  (whit-,  spectra)  or  (63),  (66), 
(89)  and  (91)  (algebraic  spectra)  one  obtains  for  remote  Interference 

rms  error  of  nulling  trac k^r _ 

mis  error  of  ccinventlor.Hl  tracker 


N  -  i _ C-l-l)'  4  fK-2)^ 

2/  2  '3N23  2 

(M-1)  \  him  +1)  -1-  ^  7  M  +  7  (-■  )  >< 

V  2  1  4  1 


(94) 


With  M  •>  1  this  reduces  to 


rms  error  of  nulling  tracker _ 9  h' _ 1 _  (“^b) 

rms  error  of  conventional  tracker  V  8  ! 

In  the  limiting  case.s  ■'f  ambient  nuise  Poninat..  :  and  inurference 
dominated  envl  ronmcnc.s  one  obtains  fJnil' 

(  I/' 2  ^  ■  ■  \/% 

rms  error  ov  nulling  tracker  _ 

rms  error  of  con  -nt tonal  tiuck-r 


(9b) 


IV 


/gM  Y  for 


l/M  -> 


.3  b  3  (N) 

2  1  9  1 


It  is  interesting  that  the  .superiority  (if  any)  of  the  njlllng  tracker 
over  the  conventional  tracker  is  tnv  asured  by  .M  ^  y  ,  the  parameter  whicii 
also  measures  directly  the  extent  lo  which  the  erwJionment  i.s  Interfer¬ 
ence  dominated.  Only  in  a  strongly  Interfi-n.ncc  dominated  environment  i.s 
there  a  significant  reduction  in  nas  i:rrot  du«.  to  nulling.  The  .slight 
inferiority  of  the  nulling  tr.acker  in  an  ambieiM  nol.se  dominated  environ¬ 


ment  is,  of  course,  due  Lo  the  f.ubop  L  ima  1  n.iture  o.f  the  primitive  inter¬ 
ference  nulling  scheine  employed  In  the  postulated  instrumentation. 


is  a  common  occurrence-.  If  this  Is  not  the  ease  one  is  interested  primarily 
In  eliminating  the  syst  vatic  error.  Since  the  latter  depends  only  on  in¬ 
terference  bearing  and  such  average  properties  a.s  interference  power  (but 
not  on  the  detailed  interference  time  function)  one  can  envision  much  .simpler 
schemes  than  nulling  to  achl.ve  the  necessary  correction. 

Even  when  the  environment  is  interference  dominated  one  m.'.y  be  con¬ 
cerned  primarily  witii  the  simplot  proble.m  of  systematic  error  correction. 
According  to  2)  the  largest  error  component  is  likely  to  be  systematic 
and  failure  to  track  may  well  be  du.-  largely  to  this  component.  Elimination 
of  the  systematic  error  alone  may  then  result  in  satisfactory  operation, 
even  though  random  error  reduction  through  nulling  or  related  schemes  might 
achieve  further  improvements  in  accuracy.  in  view  of  the  inherent  complexity 
of  all  schemes  designed  to  eliminate  time  functions  of  interference  or  other 
spatially  coherent  coroponent.s  of  the  noise  field,  it  appears  desirable  to 
study  further  whether  simple  modifications  of  the  convention;-.!  tracker, 
utilizing  only  statistical  Infoixintion  about  tlie  interference,  may  not 
achieve  satisfactory  performance. 

The  b.-.sic  difference  between  sy.stemai  Ic  and  random  error  compensnti  a 
appears  in  a  somewhat  different  guise  in  the  detection  problem.  The  av.  : - 
age  bearing  response  of  a  conventional  director  t)perating  in  the  pre.sence 
of  Interference  might  assume  the  form  sketched  in  Figure  27.  The  proximity 


'  e 

l.-'yct  Interference 

Figure  27 


.\  -  .‘.ii 


of  the  Interferenct  imposes  a  sloping  background  on  the  target  peak.  De¬ 
pending  on  the  method  chosen  to  display  the  Information  this  may  or  may  not 
affect  the  ability  of  an  operator  to  detect  the  target.  Ignoring  such  dis¬ 
play  problems,  however,  a  fundamental  difficulty  arises  only  when  the  rms 
fluctuation  becomes  comparable  with  the  target  peak,  i.e.,  when  the  "on 
target"  signal-to-nolse  ratio  at  the  detector  output  falls  to  the  neighbor¬ 
hood  of  unity.  If  the  interference-free  output  signal-to-noise  ratio  is  well 
in  excess  of  unity,  then  the  introduction  of  an  Interference  separated  from 
the  target  by  more  than  a  beamwldth  can  achieve  reduction  of  the  output 
signal-to-nolse  ratio  to  a  level  of  unity  only  if  the  interference  is  strong 
enough  to  dominate  the  environment.  In  such  cases  nulling  can  clearly  be 
beneficial.  Otherwise  the  problem  is  primarily  one  cf  display,  and  while 
this  counterpart  of  the  systematic  error  problem  Is  clearly  resolved  as  a 
by-product  of  any  nulling  procedure,  one  should  recognize  that  simpler 
techniques  may  be  available  to  compensate  for  the  average  effect  of  the  in¬ 
terference  . 
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Sunmary 

A  straightforward  extension  of  the  methods  of  Report  No.  17  is  employed 
to  obtain  general  expressions  for  the  index  of  performance  of  a  conventional 

detector  when  multiple  independent  point  source  interferences  are  present. 

-oj  It 
0 

As  in  Report  No.  17,  autocorrelation  functions  of  the  font:  c  are 

assumed  for  the  nolst  and  interference  processes,  and  the  index  of  performance 
is  derived  for  cases  where  all  interferences  differ  at  least  several  degrees 
in  bearing  from  the  target. 

An  interference  distributed  over  an  arc  of  several  degrees  is  represented 
as  a  limiting  case  of  a  large  number  of  closely  spaced  point  interferences. 
Computer  calculations  relevant  to  detemilnation  of  the  index  of  performance  In 
cases  of  multiple  or  distributed  Interferences  arc  presented. 

If  the  total  interference  power  substantially  exceeds  the  ambient  noise 
power,  the  magnitude  of  the  total  interference  power  dtctrralncs  the  order  of 
magnitude  of  the  index  of  performance.  If  a  fixed  value  of  total  interference 
power  1  is  considered,  the  case  of  a  single  point  Interference  of  power  I  may 
be  compared  with  the  cases  of  distributed  or  multiple  interferences  of  total 
power  I.  Assuming  the  average  bearings  of  interferences  to  be  comparable,  one 
finds  that  the  index  of  performance  is  somewhat  higher  for  multiple  or  distri¬ 
buted  Interferences  than  for  a  single  point  source  of  interference.  The  factor 
by  which  the  index  of  performance  with  multiple  or  distributed  interferences 
may  exceed  that  for  one  point  interference  is  found  to  have  a  hypothetical 
maximum  of  approximately  >^2M/3  ,  where  M  is  the  number  of  hydrophones.  The 
computed  results  indicate,  however,  that  the  improvement  factor  in  most  real¬ 
istic  situations  is  no  greater  than  2  or  3  for  multiple  point  source  in¬ 
terferences  and  between  1  and  2  for  a  distributed  interference. 

Finally  the  report  points  out  rhe  omission  of  a  term  in  the  basic  ex¬ 
pression  in  Report  No.  17  for  the  variance  of  the  output  of  a  standard  deteetor 
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in  the  presence  of  ambient  noise  and  a  point  source  interference.  The  term 
in  question  involves  the  product  IN  of  interference  power  and  ambient  noise 
power.  Calculations  are  performed  which  indicate  'hat  the  error  Introduced  in 
the  output  signal-to-noise  ratio  (index  of  performance)  by  neglecting  this 
term  is  at  most  a  lew  percent  for  the  values  of  N/I  and  number  of  hydro¬ 
phones  M  considered  in  Report  No.  17.  The  effect  of  the  term  is  found  to 
decrease  as  either  M  or  T/N  increases. 


I.  Introduction 


.COi  >aa  s 


The  present  report  is  essentially  an  amendment  and  extenaion  to  the  por¬ 
tion  of  Report  No.  17  which  treats  the  conventional  detector.  The  expression 
in  that  report  for  the  variance  of  the  detector  output  is  modified,  and  the 
case  ot  a  single  interference  is  extended  to  consider  any  number  of  interfer¬ 
ences  . 

This  report  considers  only  a  linear  array  of  M  equally  spaced  omni¬ 
directional  hydrophones.  Processing  consists  of  summing  the  outputs  of  the 
hydrophones,  squaring  the  sum,  and  lov-pass  ilitoring  the  square.  The  target 
signal,  the  ambient  noise,  and  the  one  or  many  interferences  are  taken  to  be 
mutually  independent  random  processes. 

The  ambient  noise  components  from  different  hydrophones  are  assumed  to  be 
independent.  We  shall  also  assume  signal  and  interference  wavefronts  to  be 
plane.  Consequently,  at  any  given  time  the  signal  and  interference  components 
of  the  M  hydrophone  outputs  represent  samples  of  these  random  processes  at  M 
different  points  in  time.  If  the  array  is  steered  on  target,  the  M  signal 
components  are  identical  at  any  time.  The  signal,  ambient  noise,  and  inter¬ 
ference  components  of  voltage  at  each  hydrophone  are  all  assumed  to  be  zero- 
mean  Gaussian  variables.  The  average  power  of  auy  of  these  distinct  processes 
is  assumed  not  to  vary  from  hydrophone  to  hydrophone-. 

The  nomenclature  of  the  present  report  is  adopted  from  Report  No.  '.7,  and 
where  specific  assumptions  are  made  regarding  magnitude  or  functional  form  of 
variables,  the  assumptions  arc  those  of  that  report.  The  ca.se  of  an  inter¬ 
ference  very  close  in  bearing  to  the  target  Is  not  discussed  in  detail. 

II.  Output  Variance  with  Ambient  Noise  and  a  Single  Interference  Present 

2 

It  is  desired  to  calculate  the  output  variance  D  (output)  for  the  sit¬ 
uation  depicted  in  Figure  1.  The  variables  i . (t)  and  n^(t)  represent  the 
Interference  and  noise  components,  respectively,  of  the  jrh  hydrophone  output. 
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This  report  follows  Che  procedure  established  in  Report  No.  3,  section  IV. 

Under  the  assumption  that  the  array  is  steered  on  target  but  that  sig¬ 
nal  power  is  negligible  compared  with  ambient  noise  or  interference  power, 
the  autocorrelation  function  of  the  output  is 


R^(t)-E  E 


i.  (t)  +  n.  (t)  E  i,(t)  +  n(t)  X 
1  h  “h  j  j  J 


E  ij^  (t-fi)  +  n^(t+T) 
k>l 


E  i.  (t+x)  +  n,  (t+x)  )  = 

i-1  ^  ^  I 


M  M  M  M 

"i  jii  kii  t'l '  ^  *  “^"’1  h"’  *  "j"’]  t*’'*"’*  * 


+  n,(t+,)]> 

Since  all  variables  are  Gaussian,  the  products  may  be  grouped  as  follows: 


M  M  M  M  I 

■  h^i  j'l  I  / "" 

[ik(t+x)  +  nj^(t+x)j  [^i^(t+x)  +  n^(t+x)]\  + 

<^|-lh(t)+nh(t)]  |-ij^(tf,)+i,,^(t+x)]\  + 


</ [ih(t)+nh(t)]  [i,(tfxHn,(t+T)]\  /|-i.(tHn^(t)]  [l^(t+xHn^(t4.x)]^  (2) 


Fig.  1 
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The  last  two  of  the  three  expressions  in  brackets  above  are  identical  ex¬ 
cept  for  Indexinp.  Since  noise  and  interference  are  independent,  the  auto¬ 
correlation  function  may  be  written 


M 


a  (t)  =  1  I  L  L 

^  h“l  j  =  l  k»l  8,=i 


ih(t)  ij(t))  +  (n^^Ct)  nj(t)^^]  [(ij^Ct+T)  i^(t+T)\ 


+  <^n^(t+T)  nj^(t+T)j>]  +  2j^/l^(t)  ij^(t+T/>  +  <^nj^(t)  n^^Ct+T)^] 

[<l.(t)  l^(t+T)^  +  <^n^(c)  nj(t+T);>j| 


(3) 


One  may  define  Che  following  normalized  correlation  functions  with  N  the 
average  noise  power  at  any  hydrophone  and  1  the  average  interference 
power • 


1  qJ^(T)  =  E  [  ih(t)  i^(  +t)]  ,  N  qJI^(T)  =  E  [n^^u)  ( t+T )  j 

In  terms  of  the  q  function,  the  autocorrelation  functioD«is 
M  M  M  M  I  ^ 

R.(1)  =  r  f  r  E  <1  [l  qMO)  +  M  q"  (0)J  [l  +  N  q|^^  (0) 


^  h=l  .i  =  l  k=l  1=1  ^ 


+  2  [l  qjj^d)  +  N  qJ^(T)j  [l  q]^(r)  +  0  q"j^(T)j  , 


The  q  terms  with  zero  argument  represent  DC  power.  Define 
as  [e^(t)  -  (DC  teims)j  ; 


M  H  M  M  j  „  , 

-t  =  ^ 

1-  A  INqJkCT)  q^,(T)^ 


> 


(M 


(5) 


(6) 


From  Eqs.  (19)  and  (li0)  of  Report  No.  3,  we  obtain  the  following  for¬ 
mula  for  the  variance  of  the  output- 

D^z)  =  I  r  r;(t)  dx  CONFiOfcN  \  & AL, 
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Now  let  and  p^(t)  represent  the  normalized  autocorrelation  function 

of  the  ambient  noise  and  intc  ranee  components  of  hydrophone  output,  re¬ 
spectively.  Since  the  noise  components  of  different  hydrophones  are  assumed 
to  be  independent,  the  following  relation  results: 


q,  .  (t)  >=  f  p  (t) 
h]  hj  n 


1  h  “  1 
0  h  j 


At  a  given  instant,  the  interference  components  at  hydrophones  h  and  j 
respectively  represent  samples  of  the  same  random  process  taken  at  an  interval 


of  seconds  apart,  where  is  the  delay  from  hydrophone  h  to  hydro¬ 

phone  k  .  It  is  therefore  true  that 


P,(T.,  +  t) 


By  Eq.  (6)  through  Eq .  (9), 


2  M  M  M  M 

D  (z)  =  ~  1  ;;  I  i:  ;”p.  (t  +  t)  d  (t  -i-  i)  dr 

h=l  j  =  l  k=l  ^ 


2N^  2  ™  2  4IN  ^  ® 

+  —  fl  P„  (T)  dr  +  —  I  I  (Tj  f  r)  M  p^(t)  dr 

h=l  j=] 


The  term  above  containing  the  product  IN  docs  not  appear  in  Eq .  (2) 

of  Report  No.  17  For  situations  where  the  ratio  I/N  is  substantially 

greater  than  unity,  this  term  does  not  significantly  affect  the  value  of 
2 

D  (z)  .  Even  with  I/N  near  or  less  than  unity,  the  term  has  little  effect 
for  sufficiently  large  M  .  Some  calculations  appear  in  Table  1  of  Sec¬ 
tion  III  which  indicate  the  relative  magnitude  of  terms  in  the  exj)ression 
2 

for  D  (z)  when  exponential  p  functions  arc  assumed. 

III.  De tector  Pe r f o rmjinco_  with  Amb lent  Noise  Two  Jn t c r f e rentes 

If  two  uncorrelated  interferences  of  power  Ij^  and  I^  are  present  at 

2 

different  bearings,  one  may  infer  the  expression  for  D  (z)  from  Eq .  (10). 
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Ambient  noise  power  Is  assumed  small  in  comparison  with  both  and  . 

„  21^  M  M  M  M  ^  . 

D  (z)  -  ~  1  E  E  E  ;*p  +  t)  p  (tJ  +  t)  di 

h-1  J-1  k-1  i-1  ^  ^ 


^  ^  "  A,  (,<P  ,)  p,  „  p, .  ol  (.)  d, 


T  u  1  ,  1  n  ,  -“2  "hk 

h»l  j=l  k^l  £=1 


2  '  ji 


<.1,  I.  M  M  M  K  ,«>  ,  (1)  ^  ^  /  (2)  ^  , 

i-  E  E  E  E  (^hk  <^2 

h=l  J-1  k-1  (l=i 


41  N  MM 

+  -T~  “  ^1^1  (^h1  +  «Pn 

h-1  j-1  J 


41  N  M  M 

dT  +  -TT^  E  E 

h=l  j“l 


^00^2  '^'hj 


In  general  terms,  p  +  f)  refers  to  the  normalized  autocorreia- 

V  hj 

cion  of  I'.he  pth  interference. 


The  procedure  of  Report  No.  17  will  now  be  followed  to  obtain  specific 

results  for  the  figure  of  merit  (D^.  following  normalized 

D  (output) 

autocorrelation  functions  are  assumed: 

-w  |v| 

P,(r)  •-  e  P,(t)  -  e  p  (r)  =  e  (12) 

X  4.  T1 

A.s  explained  in  Report  No.  37,  for  the  cases  of  interest,  the  frequencies 
(i)^,  are  comparable.  Because  the  array  is  linear  and  the  Ipiterfer- 

ence  wavefronts  plane, 

,«>  -  (b-k)  .  (h-k)  (13)1 


-Assuming  the  delays  and  to  be  positive  numbers  proves  conven¬ 

ient  in  writing  many  subsequent  expressions,  even  though  Fq .  (13)  actually 
permits  negative  I's  .  In  all  expressions  derived,  no  error  is  incurred  by 
assuming  the  r's  positive,  even  if  Eq .  (13)  indicates  that  one  or  both  are 
negative  for  the  geometry  and  indexing  of  a  particular  situation. 


For  convenience,  define  the  coefficients 


C.  .  .  =  P  dT 

h-j ,  k-i  u  hj  '  V  ki 


The  fact  that  autocorrelation  functione  are  even  functions  leads  to  the 
results 


(14) 


_uu  ^  ^uu  ^  ^uu 
rs  sr  -r,  -? 


i\v  vu  uv 

C  =  C  =  C 
rs  sr  -r,  -s 


(15) 


Confederation  of  Fq .  (7)  through  Tq .  (10)  In  Repoit  No.  17  and  the  above 
Eq.  (15)  leads  to  the  result 


H  H  M  M  „  M-1  I!-l 

1  E  E  E  C^''  +  2M  E  (M-s)  +  2M  E  (M-r) 

h=lj  =  lk=17=l  s=l 


r=l 


U 


M-1  M-1 

i  nv  UV  ' 

+21  E  (M-r)  (M-6)  I  C  +  C 

1  ,  I  rs  -rs  J 

r=l  j=l 


By  analogy  with  j.'] .  (14),  define  the  coefficients 

=  /"p  (f^V^  +  t)  p  (t)  dx 

h-j  -®  u  hj  n 


(16) 


(17) 


Using  Che  fact  that  d'’’  -  . 

r  -r 


M  M 


one  mav  show  that 


M-1 


E  E  D“  =  M  d“  +  2  E  (M-r) 

I  1  j  1  n"j  o  ,  r 

li«=l  j  =  l  r=l 

The  coefficients  are  calculated  from  Eq .  (12)  through  Eq .  (i'l).  Details 
of  the  cal  culation  appe  '.r  in  Appendix  A. 

I 

' -w  Ix+rx  I  -0)  Ix+st 

‘  I  .  '  1  1  '  >T  '  Tl 


(X8) 


rs 


u  u' 


luA  u)A  coA  u)A 

-u  -V  -u  -V 

e  +  c  o _ -e 


Ul  +  (J 
U  V 


U)  -  (J 
V  u 


(19) 


where,  by  definition, 


A  =  r  T  -  8  T 
'  u  V 


H-d 


(20) 
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If  u  ”  u  ,  the  coefficient  takes  a  simpler  form, 
u  V 


—^(1  +  (1)  |rT  -  ST  |)€ 


01  rt  -BT 
-  u  u  V 


(u)  “  01  ) 

U  V 


From  Eq.  (lA)  of  Report  No.  17, 


,  r-9  0)  T 

1  M  j  I  I  ^  "  '  u  u 

C  “ - (1  +  r-8  01  T  )e 

rs  01  u  u 

u 


By  Betting  equal  to  zero  in  Eq .  (19)  one  may  infer  the  value  of  the 

D  coefficients,  defined  in  Eq .  (17)  . 

01  I  ri  I  01  I  rt  I  01  1  ri  I  oi  I  ri 

-u'  u.  -o  u  u'  -0  u 

^u  e  +  e  .  e  -  e 

p  „ - ^ - 

r  01+0)  01-01 

u  o  o  u 


If  0'^  =  01^  ,  the  result  simplifies. 


,  oj  rx  ! 

r,u  1  /  ,  .  I  I  \  -  u  u 

D  “  -  (1  +  01  rx  )e 

r  0)  u  u ' 

u 


The  following  general  expression  for  output  variance  in  terms  of 
the  C  and  D  coefficients  represents  a  combination  of  Eqs.  (11),  (14), 


(16),  (IV),  and  (18) 


2  f  M-1  M-1  M-1  - 

D  (z)  -  -tr^{  M  C"  +  4M  1  (M-b)  C  +  2  1  I  (M-r)  (M-s)  C  +  C 

^  s-1  r-1  s=-l  L 

2l/  I  2  22  22  ■  22 

+  C  +  AM  I  (M-s)  +  2  E  E  (M-r)  (M-s)  C '  +  C 

T  )  00  T  08  1  1 

8*1  r*l  S"1  L  J 


2N^  „2  «  2  ,  ,  j  ^  ^^1^2  ]„2  ^12  ^  ,,,  v  4. 

+  — —  M  /  p  (x)  dx  +  — r — <  M  C  +  2M  E  (M-q)  >  C  +  C 

T  n  T  )  oo  ^  oq  qo 

a“i  0 


''  ^  .  [  12  12  ll  *^1^  <M^D^  +  2M  E  (M-r)  i 


4I2N  ,  ,  M-1  2 

+  —  <  M  D  +  2  E  (M-r)  D 

^  °  r-l  ' 
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For  the  assumed  exponential  autocorrelation  functions,  by  Eqs.  (19)  and 
(22)  through  (25) . 


4  4H  £  ^  (1  4  ,  , 

s»l  “^l 


H-1  M-1 

2  r  E 

r=»l  s»l  '^1 


(|+|r-s|w^Tj^)e  ^^1  + 


(l+|r+s|uj^T^)e 


■(r+s)^^!^ 


21, 


I  {terns  like  those  for  1.,  olth  replsced  by  + 
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H-l 
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1  2M  ‘  ‘■ 

~  +  2K  £  (M-r) 

“1  ‘^o  r=l 


“l‘'“o 


%"‘^1 


4I2K 


^  -(terms  like  those  fur  ,  with  replacec 


by  U2T2 


In  the  above  expression,  a  and  £  are  defined 


-r-  1  rv  ^ 

tK  '.  (26)1 

t  ' 
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|rT^  -  ST^I 


T.  -  rt^  +  sT^ 


(27) 


If  neither  Interference  is  near  the  signal  in  bearing,  the  products 

and  “2^2  much  greater  than  unity,  and  most  of  the  terms  in 

Eq.  (26)  may  be  neglected. 

2  2 

D^(z)  =  ~ 


2  1  ' 

-  M  +  ^  K| 

j-  II  m' 

2  1  V 

-  M  +  ^  Ml 

“1 

“1 

3  ; 

T 

"2 

“2 

3  J 

2N 


2  1 


T  ((j)j^+a)2) 


2 

2M  + 


2  3  2  1 

I  M  -  M  +  ^  M 


12 


T(a),+u.  ) 
1  o 


2M 


4I,N 

0 

+  - = - 

T(u..,+u  ) 

^  0 

2M‘' 

. 

(cOiT^,  W2T2  »  1) 


(28) 


The  "interaction  factor"  X^2  measures  the  effect  on  output  variance  of 
the  Intermodulation  of  the  two  interference  processes  which  results  from 
the  squaring  operation  in  the  otandard  detector.  If  both  interferences  arc 
remote  in  bearing  from  the  target,  the  factor  X22  has  a  maximum  value  of 
2  .  This  value  is  reached  when  both  interferences  have  the  same  bearing. 
The  factor  is  discussed  in  detail  at  the  end  of  this  section. 

The  effect  of  a  signal  having  power  S  when  the  array  is  steered  on 
target  is  demonstrated  in  Sq .  (33)  of  Report  No.  3.  That  effect  is 

A  (DC  output)  -  M^S  (29) 


From  Eqs.  (28)  and  (29),  one  may  calculate  the  following  figure  of  merit 


For  the  assumed  values 
feet,  calculations  show  mat 
'ngs  greater  than  about  10° 


u)  >»■  2r  x  5000  and  hydrophone  spacing  d  =  2 
the  term  e  ■*  is  less  than  .1  for  bear- 
relative  to  the  targe 


■L- 
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For  slnplicity,  assui..i.  a;,  =  u-  -  ■  then  Fq .  indicates  a  slit- 

±  I  o 

pliflcatlon  in  the  fom  of  the  terns  involving  A  . 


H-1  M-1 

2  1  Z  (M-r)  (M-s)  —  (1+0)  A)  e 
11  0)  o 

r=l  S“1  o 


^  1  't  -1  O  1 

°  .  —L-!  2. 

20,  3  '•  '^12 

01-  J 


Assune  that  =  ki^  ,  where  k  ^  1  .  The  double  suenation  nay  be 
written 


M- 1  '!“1  -0)  T  1  r-ks  1 

2  Z  L  (M-v)  (H-s)  —  (l+o)  T  ir-ke!)  e  °  (; 

r=l  s=l  “^o  °  ^ 

-W  T 

o  1  ^■*0  ^  2 

Consistent  x^ith  the  assunption  that  e  ?  e  s  0  5  one  nny 

neglect  terr.'j  in  the  sumr.ation  except  those  for  which  |  r-ks|  <'^1  .  For 

concreteness,  one  nay  require  jr“ks|  <  1/2  as  the  condition  for  a  tern 

to  be  significant.^  For  a  particular  integer  s'  ,  suppose  that  the  in¬ 
teger  r'  satisfies  the  relation  r'-ks'  a  1/2  .  Then  (r'-l)  -  ks'  = 

-  1/2  .  Ir  this  case,  for  each  s  there  are  two  values  of  r  such  that 
jr-ks|  _<  1/2  .  Except  for  a  discrete  set  of  values  nf  k  ,  no  r  will 
be  found  for  which  r-ks  =  1/2  exactly.  However,  as  long  as  k  <  (M-1)  , 
sore  one  value  cf  r  can  be  found  for  sufficiently  snail  values  of  s 
such  that  |r-ks  !  <  1/2  .  Fence,  in  r.eneral,  a  reasonable  assunption  is 
tliat  for  each  value  of  s  ,  only  one  value  of  r  need  be  considered. 

Note  that  by  definition  i  =  (d  sin  0)/c  ,  where  d  is  hydrophone 

spacing  and  c  is  sound  velocity.  Hence,  k  -■  /  i  =  sin  0^  /  sin  0^  . 

-u  T 

Since  the  basic  assunption  e  =  0  is  poor  for  angles  less  than  about 


1  I 

The  assunption  being  nad>;  is  that  j  1+u)^Tj^(1/2)J  c  °  is 

<:  1  ,  on  the  order  of  .01  ,  for  Instance.  ".-is  condition  is  true  for 

2  3.5  •  for  u,^  -  2ti  X  5000,  d  •=  2  ft.,  the  inequality  holds  for  any 
relative  bearing  greater  than  about  15°  . 


H-J  1 
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15  ,  the  inininun  value  of  sin  6^  for  which  the  analysis  is  valid  is 

about  .25  .  Hence,  k  is  not  nore  than  1/.25  ■=  4  . 

In  the  expression  (36),  s  values  greater  tha.n  (M-l)/k  find  no 
corresponding  r  satisfying  the  requirement  jr-ks|  <  1/2  .  Hence,  the 
naxirr.ur.  value  of  s  which  is  significant  is  the  Inte.ger  nearest  (M-1) /k 
(or  M/k  for  simplicity) .  The  value  of  r  corresponding  to  each  s 
is  the  integer  nearest  ks  .  The  expression  (36)  is  then  roughly 

[m/kj  -u)  T 

2  I  (M-ks)  (M-a)  —  |r-ks|)  e 

,  U)  '  o  1  ' 

S  =  1  0 

where  |  M/kj  denotes  the  largest  integer  smaller  than  M/k  .  Investi- 

g.ation  reveals  that  for  a  given  s  ,  the  value  of  i r-ks |  is  the  dif¬ 

ference  between  sd  and  Che  nearest  Integer,  where  d  is  the  fraction¬ 
al  part  of  k  .  In  general,  the  .above  expression  cannot  be  significant¬ 
ly  simplified.  The  expres.sion  has  local  maxir.ia  with  respect  to  k  when 
k  is  a  multiple  of  1/2  ;  local  minima  occur  approximately  where  the 
fractional  part  of  k  is  1/4  or  3/4  . 

A  rough  estimate  of  the  summation  nay  be  obtained  for  the  cases 
where  k  is  an  integer.  For  all  values  of  s  yielding  significant 
terms  in  the  sumnation,  |r-ks|  =•  0  .  Hence,  by  (37)  the  summation  re¬ 
duces  to 


|r-ks! 


L  m/k 
2  I 
s=i 


(M-ks)  (M-s)  --  i  —  ;  ™''’'(M-ks) 

0'  (J  i 
o  o 


(M-s)  ds  = 


2  m/k  2 
^Oo  1  L 


-  M(k+1) 


s 


K^s- 


(k+1) 


3 

s 


n/k 

1 


For  the  product  uj  d  on  the  order  of  4it  x  5000  . 
'  o 
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(37) 
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2  il_  _  „2  _  fl„(k+ll  ,  M  ^  ^  _  _k 
“  k  ,,2  '  2  ^  „3  3 


2  ^3  ,1  k+1  ,  1 

—  M  (-  -  — -  +  — 

0  ^  2k  3k 


m21  _  2 

v2^  *  ^  7 

^  -“ok 


(k  -  H  (k+1)  +  3) 


^  m:  k  _  1 

a  ,  2  '■2 
0  k 


For  large  M  ,  the  above  expression,  where  powers  of  M  less  than  the 

third  power  are  ignored,  should  be  reasonably  accurate-  Pcferring  to 

2  3 

(B3) 3  and  assuming  M  <<M  ,  one  obtains  an  approximation  for  • 


-k/2  ■  1/6  _  3k  -  1 

k^ 


=  integer) 


It  appears  Impossible  to  obtain  a  simple  analytical  expression  for 
the  factor  Xj^2  terms  of  the  bearings  0^  and  02  and  the  para¬ 
meters  of  the  system.  For  this  reason,  extensive  computer  calculations 
have  been  performed.  For  the  calculations,  it  is  assumed  that  both 
processes  have  the  same  bandwidth.  The  precise  form  of  the  expression 
for  X^2  >  derived  from  Eq .  (35)  Is' 

N- 1  M-1  -m  I  ri.. -ST2  | 

X  =  -nr----r--  f  I  (M-r)  (M-2)  (1+m  jrt  -sT J)  e  ° 

2M  -3M  +M  r=l  s=l  o  i  / 

M-1  M-1  oj  d 

=  — 2 - 3 —  2  T.  (M-r)  (M-s)  (H— ^  Ir  sin  0  -s  sin  02l)X 

2M'^-3ir+M  r-1  s=l  ^ 

ii)  d 

O  I  ,  ^  r  .  I 


|r  sin  Oj^  -  s  sin  O2  1 


where  M  is  the  number  of  hydrophones,  d  is  hydrophone  spacing,  and 

c  sound  velocity.  For  calculation,  m  v;as  taken  to  be  2^  x  5000, 

o 

d  two  feet,  c  5000  feet  per  second.  Results  were  obtained  for  M  =  40 


and  M  =  20. 


jV.'i )  iri: 


In  Figure  2  curves  of  versus  0  ,  the  arithmetic  mean  of  the 

bearings  of  two  point  source  Interferences  relative  to  the  target,  are 
plotted  for  different  constant  values  of  a0  ,  the  difference  in  bearing 

between  the  interferences;  for  these  curves  M  is  taken  to  te  iO  hy¬ 
drophones.  In  Figure  3  curves  of  versus  breadth  ^0  ,  for  fixed 

values  of  center  angle  0  ,  are  plotted  for  both  M  =  40  and  M  "  20 
to  permit  comparison  of  the  results  for  different  values  of  M  . 

In  Figure  4  contours  of  constant  delay-time  ratios  are  plotted 
on  the  0^  -  0^  plane. 

A  striking  feature  of  the  curves  (Figure  2)  is  that 

for  a  given  value  of  0  in  the  range  of  about  25°  to  65°  ,  X^2  fs 

nearly  independent  of  Che  breadth  (angular  separation)  A0  for  values 
of  AO  between  about  5°  and  40°  (.1  radian  to  .7  radian).  Exceptions 

to  this  statement  do  occur  in  certal  •  ranges  where  one  or  more  curves  have 

sharp  relative  minima  or  maxima  with  respe:t  to  0  .  The  limited  calcu¬ 
lations  for  AO  on  the  order  of  one  radian  indicate  t  at  X^2  does  dip 
nearly  to  zero  for  at  least  one  range  of  0  • 

The  curves  may  actually  be  somewhat  more  Irregular  than  the  plots 
Indicate,  since  the  points  picked  for  calculation  did  not  include,  all  the 
combinations  of  angles  yielding  delay-time  ratios  which  correspond  to 
relative  maxima  or  t.inima  of  X^2  •  It  seems  reasonable  to  conclude,  how¬ 
ever,  that  unless  0  is  fairly  near  90°  or  AG  is  on  the  order  of  a 
radian  or  greater,  detector  performance  is  relatively  insensitive  to  the 
value  of  AO  .  It  must  be  noted  that  for  M  set  equal  to  20  ,  the 
curves  separate  appreciably  at  a  smaller  value  of  0  than  for  M  equal 
to  40  .  Evidently  a  larger  number  of  hydrophones  Increases  the  range 
of  0  over  which  performance  is  insensitive  to  AO  . 

The  comparative  flatness  of  the  curves  of  X^2  versus  AO  (Figure  3) 
illustrates  the  above  remarks. 
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IV.  /^jnendments  i.o  Re  port  ’-'uinber  17 
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On  the  basis  of  the  derivations  in  the  previous  section,  it  will 
new  be  indicated  where  in  i.eport  No.  17  terms  need  to  be  added  to  ex¬ 
pressions.  For  direct  comparison  wdth  the  result^  of  the  earlier  re¬ 
port,  it  must  be  assumed  that  the  interference  "cutoff"  frequency 
is  the  same  as  .  By  Eqs.  (il),  (17),  (18).,  (2A),  and  (25),  Eq .  (16) 

of  Report  No.  17  Is  amended  as  follows; 

A(DC  output)  ^ 

D(output) 

/t„i  f,;g 

[  /  _ O - - 

2  r  .  .  f  ."-1  >  T  M-1  .  -2sm  - 


-BW  T  M-1 

00.2  .  , 


N  +I  (  1+--  1  (M-s)  (l+S(o  T  )e  +  --q-  L  (.M-s)"  .l-!(l+2suj  t  )e 


M  , 
(  .s=l 


:.’-l  r-1 


—  1  1  (M-r)  (M-s)  l+(r-s)o)  t  'e 

r=2  0=1  U.  ° 


M  H=1 

•1  -(r-s)ii  1  r 

0  0, 

'  e  +  . 


-2su)  T  I 

0  o 


l+(r-t-s)u)  T  e 
o  o 


■(i4-s)i..  r  \ 
o  o  I 


/  ,  M-1 

+  IN  (  2  +  ^  I  (M-r)  (1+ru  t  )e 
M  ,  0  0 

I  r-1 


-r„  ,  1 

o  o  \ 


Accordingly,  Eq .  (17)  of  Report  No.  17  is  amended  to 


i(DC  output)  ^  _ 

B(output)  V  2  j - 

/  2  2 

.  /  N  +1 


I  1  +  ^  >-  (M-s) 

M  s=l  -I 


+  2IN 


Eq .  (18)  of  Report  No.  17  becomes 


D(outpuL) 


/  2  I  1 

/  N  +  Y  (2M  +  ^)  +  2IN 

In  terms  of  the  factors  dc  .fined  in  Eq .  (32),  tlie  above  equation  reads 


A  (DC  output) 
D(output) 


\f  N^+A1^+21N 


Fq .  (A3)  reveals  that  the  term  2IN  is  relatively  unimportant  unless 
M  and  I  are  comparable  and  M  is  not  large.  To  determine  precisely 
the  effect  of  this  term,  one  may  assume  I  =  N  .and,  using  Eq .  (AA)  , 
calculate  a  correction  factor  ns  a  fund  lor.  of  M  which  will  convert 
Eq .  (18)  of  Report  17  to  agree  with  Eq .  (A3)  above. 
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Table  1  (N  -  ") 


_M 

Correction  factor 

10 

.89 

20 

.94 

30 

.96 

40 

.97 

50 

.97 

Since  some  of  the  curves  In  Report  No.  17  take  M  as  40  and  allow 
N/I  to  vary,  a  correction  factor  as  a  function  of  N/I  may  be  useful. 

Table  2  (M  «  40) 


m 

Correction  factor 

.1 

.99 

1 

.97 

4 

.92 

V.  Detector  Performance  with  Ambient  Noise  and  Several  Interferences 

Without  further  derivation,  one  can  infer  the  expressions  for  output 
variance  and  figure  of  merit  for  k  Interferences  from  Eqs.  (25),  (28), 
and  (30).  Using  the  C  and  D  coefficients  defined  In  Eqs.  (14)  and 
(17) ,  the  general  expression  for  output  variance  is 
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1  ^  L  oq  qoi  L  rs  -rsj 


+  ^  E 


.2_u 


M-1 
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M  D“  +  2M  E  (M-r)  D“  ^  ^  N  <  M  (r)  dx 

r«l 


(45) 
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If  autocorrelation  functions  are  assumed  to  be  exponential  functions  of 

the  form  e  and  all  products  are  much  greater  than  unity, 

the  index  of  performance  becomes 

A (DC  output)  ^ 

D (output) 


'f 


f  /  \ 

k  I  2  k 

k  II 

- 

j  1  2L  1> 

r  **  4.  r 

U  V 

2+1  |m-H 

1^^  1  y 

Z  - +  E 

U-1  *^u  U“1 

1m  i  uv 

k 

+  2N  E  ■ 
u=l 


I 


u 

is(aio+a.^) 


(46) 


where  is  the  same  function  as  X^2  Section  III,  Involving  the 

uth  and  vth  interlerences. 

It  seems  instructive  to  compare  detector  performance  in  the  presence 
of  several  independent  point  source  interferences  with  the  performances 
in  the  presence  of  a  single  Interference  yielding  the  same  average  power 
as  all  the  Independent  interferences  together.  If  I  represents  the 
total  average  interference  power  in  both  situations. 


I  = 


k 

r 

u=l 


I 

u 


(47) 


2 

In  Eq.  (46)  the  coefficient  of  the  I  ierms  is  the  factor  A  de- 

u 

fined  in  Eq.  (32),  and  the  coefficients  of  the  I  I  terms  are  like  B 
of  Eq.  (32).  B  ranges  from  2  to  2A  as  Xj^2  varies  from  0  to  2  . 

In  the  limiting  case  where  angular  displacements  between  interferences  are 
small  (or  bandwldths  are  narrow),  the  B  factors  in  Eq.  (46)  approximately 
equal  2A  .  If  all  the  roughly  equal  ,  Eq.  (46)  then  becomes 


■^ote  that  2.0  is  the  maximum  value  of  X^^.,  for  Interferences  remote 
in  bearing  from  the  target.  For  less  remote  Interferences  X.j  exceeds  2  . 
Wherever  the  assumption  e  °  ^  s  0  is  valid,  X^2  ®1^  ^  ' 
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1  ^ 
U=1 

1  ^ 

1  u*=l 

AI^  +  2rN  + 


(48) 


This  result  is  the  same  ns  that  obtained  in  Eq .  (44)  for  a  single  inter¬ 
ference  of  power  I.  it  seems  reasonable  that  if  multiple  interferences 
span  only  a  small  arc,  on  the  order  of  the  beamwidth  of  two  adjacent  hy¬ 
drophones,  the  effect  on  the  detector  is  scarcely  distinguishable  from 
that  of  a  single  point  source  Interference. 

The  opposite  limiting  case  occurs  In  the  very  improbable  event  that 
the  interferences  are  all  located  at  critical  bearings  which  make  nil  the 
approximately  zero.  In  this  situation,  the  B  factors  roughly 


The  best  performance  smallest  D(output)]  occurs  if  all  the  are 

equal  to  lu  this  event, 


This  conclusion  results  from  minimizing  the  denominator  of  Eq .  (49) 
\.lth  the  constraint  El  =  I  . 

CONFIDENTIAL 

B-21 


V  2 


D(output) 


T,^ 


MS 


Tu) 


MS 


0  f  fi  ^  9  ^  ^ 

I  ,NI  +  H  ^ 


+  2NI  +  N 


For  3  ratio  N/I  1  ,  this  result  is,  of  course,  about  the  same  as  that 
for  a  single  Intcrfun-'n'i  .  If  1  >>  N  ,  however,  the  figure  of  merit  in 
this  case  Ir.  greater  than  that  for  a  slng.^e  interference  by  a  factor  of 


Ak 


about  /  If  k  «<  10  and  M  »  40  ,  for  instance,  this  factor  is 

A■^k 


roughly  3  .  In  practice  this  performance  would  virtually  never  be 
achieved,  since  it  depends  on  a  freak  distribution  of  the  interferences 
in  space. 

A  crude  but  hopefully  more  meaningful  estimate  of  the  best  perfor¬ 
mance  for  a  fixed  total  interference  power  is  obtained  by  assuming  that 

all  the  X  in  Eq .  (46)  take  on  the  minimum  calculated  values  for 
uv  ^ 

1-^)  with  AO  in  the  range  of  about  5*^  to  40°  .  This  minl- 

6  .  From  Eq .  (46) , 


mum  value,  for  M  =  40  or  M  =  20  ,  is  about 
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MS 
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(X  =  .6,  all  u,  v) 
uv 


With  large  I/N  ,  this  figure  of  merit  is  greater  than  that  for  a  single 

For  k  -=  10  ,  this  figure  la 


interference  by  a  factor 


l+.3k 

about  1.6  ;  for  large  k  ,  it  is  about  1.8  . 
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A  crude  estimate  of  average  performance  is  obtained  by  letting  all 
the  in  Eq .  (46)  take  on  the  average  calculated  value  of 

A0  In  the  range  of  5°  to  40°  .  The  average  is  close  to  1. 


A(DC  output) 
D(output) 


fi  *  f) 


+  2NI  +  N 


(X  -■  1.  ,  all  u,  v) 
uv 


/  k 

The  improvement  lactor  here  is  ,/q - —  ,  which  is  about  1.2  for 

V  1+ .  5k 

k  =  10  ,  and  about  '/  2  for  large  k  . 


VI ■  Detector  Perform ance  with  Ambient  Noise  and  a  Dlstrlbi ited  Interference 

An  interference  distributed  continuously  over  a  finite  arc  may  be 
treated  as  a  limiting  case  of  the  multiple  interference  problem.  The 
distributed  inter ferenct:  is  represented  as  an  infinite  number  of  elemental 
point  source  interferences  spaced  an  infinitesimal  angular  distance  apart 
throughout  the  arc.  The  derivation  in  the  previous  sections  of  this  re¬ 
port  have  assumed  the  several  interference  processes  to  be  statistically 
independent,  and  the  results  to  be  derived  here  will  not  reflect  depen¬ 
dencies  -among  different  points  along  the  arc  of  the  tlistributed  interfer¬ 
ence.  Tlicse  dependencies  may  be  expected  to  degrade  pe'..  f  ormance  somewhat 
more  than  the  subsequent  results  of  ehis  section  will  indicate.  The 
assumption  of  expone.itlal  autocorrelation  functions  for  interference  and 
noise  is  inherent  in  these  deriVr.Llotis,  and  the  'cutoff  '  frequencies  for 
interfertnre  and  noise  arc  assumed  equal. 

Angular  power  density  functions  may  be  defined  ,as  follrwt  to  cor¬ 
respond  to  the  terms  of  Eq .  (46); 


1(9  )  d0  »  I 

u  u  u 


i(e  )  do  ■=  1 

v  V  V 


(53) 


The  first  two  expressions  In  the  denominator  of  Eq.  (38)  may  be  combined 
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and 


X  (e  ,  ,  e  ) 

min  max 


(9  -0  ^ 

max  min 
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max  max 

72  e  ti 

)  min 


de  x(e  ,  6  ) 

V  U  V 


(59) 


Now  Eq.  (A6)  is  approximately 


A (DC  output) 
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r  T  1  ^  x(®  .  6  ) 

'  r2  i-i.liMi4.i-  min*  max 

I  1  +  1  3  M-1  +  J  2 


(60) 


+  2NI  +  N 


The  factor  X  may  be  calculated  numerically,  but  analytic  approxima¬ 
tion  does  not  appear  feasible.  Computer  calculations  have  been  performed, 
using  the  parameter  values  ■  2ii  x  5000,  d  »  2  ft.,  and  both  M  -  40 
and  M  «  20.^  The  results  are  shown  graphically  in  Figures  5  and  6  in  the 
same  format  used  earlier  for  X^2  ■  The  curves  of  X  resemble  those  of 
Xj^2  1  except  that  the  X  curves  show  no  significant  irregularities,  and 
X  docs  not  approach  0  even  for  large  values  of  the  parameter  46  (full 
angular  spread  of  the  Interference).  It  must  be  emphasized  that  the 
factor  X  and  the  corresponding  detector  performance  are  quite  insensi¬ 
tive  to  the  breadth  of  the  interference  if  tne  center  of  the  interference 
is  closer  to  the  target  than  about  50°. 

The  results  are  not  drastically  different  lor  the  two  values  of  M  , 
the  number  of  hydrophones.  For  a  fixed  center  angle  0  ,  the  X  curves 
for  M  «  20  lie  above  those  for  M  »  40  at  small  v-alues  of  46  and  then 
approach  the  curves  for  M  «  40  at  larger  values  of  46  ,  The  difference 
in  results  for  the  different  values  . f  M  is  seen  to  be  more  pronounced 


Details  of  the  program  in  Appendix  B. 


n-25 


An 


CONFIDENTIAL 

DUNHAM  LAOORATORY.  YALE  UN1vet?SlTY 


Form  E£>12 


for  large  values  of  breadth  A6  .  The  implication  of  these  results  is 
that  an  increase  in  the  number  of  hydrophones  effects  a  small  improve¬ 
ment  in  relative  performance  as  measured  by  the  figure  of  merit.  For 
the  smaller  number  of  hydrophones,  as  compared  with  the  larger  number, 
performance  is  degraded  more  seriously  by  angularly  narrow  interferences 
than  by  broad  interferences. 

As  in  the  previous  section,  one  may  compare  performance  in  the 
presence  of  a  distributed  interference  of  total  power  I  with  that  ob¬ 
tained  for  a  point  source  interference  of  the  same  power.  For  the  case 
where  I/N  >>1  ,  comparison  of  Eqs.  (A4)  and  (60)  indicates  that  the 
'improvement  factor"  for  a  distributed  interference  is  roughly  /  2A/B  , 
where,  by  analogy  with  Eq .  (32), 


Table  3  below  displays  sample  results  for  40  hydrophones. 

Table  3  (distributed  interference) 

Center  Improvement  _ 

.Angle  0  Breadth  AO  X  B/2  A  Factor  /  2A/B 

(degrees) 

(d.-r,r..e  s) 

25 

5.7 

1.60 

22.4 

26.7  1.09 

11.4 

1.58 

22.1 

1.10 

22.9 

1.58 

22.1 

1.10 

40.0 

1.58 

22.1 

I.IG 

50 

5.7 

1.18 

16.8 

1.26 

11,4 

1.02 

14.6 

1.35 

22.9 

.95 

13.7 

1.40 

40.0 

.95 

13.7 

1.40 

75 

5 . 7 

1  .51 

?1  .2 

1.13 

11.4 

1.21 

17.2 

1.25 

22.9 

1.00 

i  ^ 

1.36 

40.0 

.97 

U  .  0 

,  1.38 
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The  "improvement  factors"  indicated  in  Table  3  do  not  reflect  a  great  dif¬ 
ference  in  performance  between  the  two  cases  of  a  point  source  interfer¬ 
ence  and  a  distributed  Interference.  Yet  a  ! ac'  jr  of  1.3  or  1.^  is 
not  entirely  trivial,  because  of  Che  fact  that  even  if  X  approached 
zero,  B  could  be  no  smaller  than  2  ,  and  hence  the  Improvement  factor 
1  nild  b  no  larger  than  A  ,  which  in  this  case  is  5.2  .  The  meaning 
;  the  hypothetical  Improvement  factor  A  becomes  clear  when  one  con- 
.siders  !>; .  ;  the  single  interference  case,  which  reads, 


»¥>■«  h 

discrepancy  between  1.5  and  wou.ld  then  have  to  be  attributed  to  the 

phone  to  phone  noise  coherence  which  exists  even  in  the  isotropic  case  with 
the  numerical  parameters  considered. 

VII.  Ccncludlnp,  Remarks 

Ai.thcugh  calculatioiij  of  multiple  or  distributed  in- 

Lerferences  may  become  quite  complicated,  the  upper  and  lower  bounds  may 
always  he  quickly  determined  if  the  total  interference  power  1  is  known. 

The  upper  and  lower  bounds  correspond  respectively  to  treating  I  as 
'fotropic  ambient  noicc  independent  from  hydrophone  to  hydrophone  and  as  a 
pi' InL  source.  The  following  inequality  must  be  satisfied. 

/  ’^'^o  MS  a(DC  output)  MS 

/  '"2"  "7==:======:  <  ^777^0“  ~ 

i/  '  -n 

/  1  I  M  j 

The  nuniericul  results  oi  this  report  indicate  that  in  r:ost  realistic  situa¬ 
tions  the  index  of  performance  will  be  much  closer  to  its  lower  bound  than 
to  its  upper  bound. 

The  major  failing  of  the  treatment  given  in  this  report  is  that  the 
results  are  inaccurate  for  interferences  ne;ar  the  target  in  bearing. 
Approximate  computations  and  the  results  of  Report  Mo.  17  suggest  that  the 
error  will  be  small  If  all  interferences  are  separated  from  the  target  by 
appreciably  more  than  the  beamwldth  of  the  array. 
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Appendix  A:  Calculation  of  C  Coefflcienls 

Initially  one  may  assume  r  grea'er  than  sr^  ,  both  r  ard 
are  positive.  Now  Eq .  (19)  may  oe  expressed 
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uuvv.uuv,  uuvv,.  V  uv, 
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-or  _  P  X 
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Lsult  generalizes  to 


Appendix  B:  Computation  of  X^2  ^ 

The  computer  program  .jnployed  to  calculate  X  and  X  embodies  a 
sub-routine  which  calculates  Xj^2  ^  function  of  6^  and  62  according 
to  Eq.  (40),  with  the  minor  exception  that  the  term  M  Is  neglected  In  the 
constant  factor. 

The  factor  X  Is  a  two-dimensional  Integral  over  a  square  In 
space  with  6^  and  6^^  running  from  6^^  to  62  •  The  main  program  por¬ 
tions  the  square  region  Into  one  hundred  small  squares  of  equal  size.  The 
value  of  X^2  center  of  each  small  square,  calculated  by  the  sub¬ 
routine,  Is  used  as  the  value  of  X^2  whole  square.  The  approx¬ 

imation  to  X  Is  then  just  a  Rlemann  sum  based  on  one  hundred  squares. 

This  method  of  approximation  yields  results  which  one  expected  to  be 
slightly  larger  than  the  correct  value  of  X  because  of  the  fact  that  at 
the  center  point  of  ten  of  the  one  hundred  squares,  ~  \ 

6  ■  0  ,  X, -  (0  ,0  )  has  a  value  of  at  least  2.0  .  For  small  differ- 

ences  (0^  -  0^^)  Xj^2  falls  fairly  rapidly  from  its  value  of  2.0  or 

greater.  Hence,  the  approximation  to  X  Is  too  large  on  these  ten  squares. 

In  fact  the  result  of  the  approximation  to  X  cannot  be  less  than 
.1  (2.0)  B  .2  .  All  of  the  calculated  results  cited  in  this  report  are 
substantially  greater  chan  .2  ,  and  it  Is  believed  that  Che  error  in 
approximating  X  Is  at  most  a  few  percent.  Spot  checks  were  made  using  a 
higher  precision  method  which  partitioned  the  0^  -  0^  space  Into  squares 
.01  radian  on  a  side.  The  results  obtained  in  this  way  were  lower  by  about 
five  percent.  Since  the  higher  precision  calculation  required  a  long  ex¬ 
ecution  tine,  only  a  few  calculations  were  made  with  it. 
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